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Stetloa 1.0 


INnODUCTZON 


Tba taeeattfai 4aslta of theraol protection cyetems for ▼chieloc operatiag 
la ataosphere and near>epaee oavlroaneatc re^airec accurate aaalyaea of keatiag 
rate and teaperature blctoriea encountered along a trajectory. For prellainary 
deelgn ealcnlationa, however, the requlreaent for accuracy auet be teapered by 
the need for apeed and versatility in coaputational tools used to deteraine such 
**tberaal enviroaaents.* * 

The NINIVES prograa (Ref* 1) over the last decade has been found to pro<- 
vide the proper balance between versatility, speed and accuracy for an aerother~ 
aal prediction tool* The advaneeaent in coaputer aided design concepts at Lang- 
ley Research Center (LaRC) in the past few years has aade it desirable to incor- 
porate the MINIVER prograa into the LaRC AVID system (Ref* 2) * The purpose of 
the AVID systea is to provide the prellainary design engineer with a useful tool 
for aulti-discipline interactions to perfora partial or coaplete vehicle syn- 
thesis* 

In order to effectively incorporate MINIVER into the AVID syatea, several 
changes to MINIVER were aade* The thermal conduction options in MINIVER were 
removed and a new Explicit Interactive Iheraal Structures (EXITS) code was de- 
veloped* Naqr upgrades to the MINIVER code were made and a new Langley version 
of MINIVER called LANNIN was crested* 

This report is divided into three volnaes. Voluae I describes the theoret- 
ical methods and subroutine functions used in LANNIN* Voluae II provides a user 


»•?« ,.id. tor LANHIN. V,l„. 1,1 4...r,b.. u. EXITS ..4, ,.4 ,.oW4.. 
iapot tnid*. 


n. 4..M.M.U.. ,r.,..„4 i, VolM. I „m.,4 ,„0, 4o.»..„„.. ,l,„, 
.»»u«»l. ..4 .44.4 M. TE. ,, 

prior 4o«uB..t.Uo. ... t.{, 1 tor tho ..Isiul EINIVER co4«, I.f. 3 
b.o.tlM 4..«H,,l.n, oot Bot. 4 tor rorotUi tlo, i;4.t... 


8«fitio& 2.0 


GENERAL UTILITX ROUTINES 


Tb« proitMi «••• tvbroatlnts of losortl ntilitp oocli •• tpooifica*” 

tloft of priatlat of iapat oatpnt data, istorpolition toelmiqMt and tho atar 
apiat of all othar oabxontiaaa. This saetioa gives a descxiptioa of these atil** 
ity xoatiaes with a dieeassioa of their ospsbilitiee. 

2.1 HABI 

MAIN perfoxBS the task of suiaagiag the rest of LANMIN. This roatiae ia 
best described by the flow diagraa ia Figure 2-1. The MAIN roatiae sets up the 
iaput. ealls the esleulatioa subroutiaes ia proper sequeueo. saves sad eoatrols 
the output, aad eoatrols the esse to esse iaput. Most of the iaput ia aeeoai- 
plished through aa array 1(700). Msay of the iaput array eleaeate are 
equivaleat to other variable aanes. A eoaplete deseriptioa of the iaput is 
givea ia Voluae II of this report. 

2.2 SDBROOTINB 1R1MP 

Sabroutiae WRINP is eslled by MAIN to priat out all of the iaput data for 
eseh separate ease. Ikis routiae priats out the pertiaeat iapat iaforastioa. 

2.3 SDBRODTINE VANOOT 

Subroutito. 7AN0UT is eslled by MAIN sod is used to priat the output ia a 
apeoified foraat. The paraaeters ia the CALL stateaeot of VANOUT are aa 
follows: 

D - (iaput) This astriz diaeasioaed (100. 13) eoatsias 13 variables as 
a fuaetioa of tiae. The 13 variables are: tiae. altitude, velocity. 

Mseh ouaber. aogle of sttsek. Reyaolds auaber/ft. heat traasfer coef- 
fieieot. recovery eathslpy. equilibriua wall teapersture. eoaveetive 
hestiag rate, heatiag load* loesl pressure, sad flow type. 

IP - (iaput) Nuaber of tiae poiots. 


<»•) «<..t.l.. .. „tu 

Mi?r’ 2 etI*J‘JI“* •» W(« 4 «) 4.t,r.U.t th> 

•IN - (1>PM) Bo <7 pout ,»bt. corrt.po.dip, to 

ni. ..»ro.tl« writ., to p.p.r o.d to . file. p..b ,, 

oitb . u,io, . 

2 .4 SOfiRODTlNE 1BLIN 

Sttbro.tlM mw i> 0 Otuitp ll.„r l.t.rpoUtio. ro.tUo oki.b u 

o.to„ivo.P tbro.,bo.t UN«N bp otb.r rootio... ,„uooi.. di.o...,o. oo.- 

cidor. tb. „o Of tbl. ro.ti« oopjetio. „tb .. .t.o.pb.ri. t.bU .blob 

-. 0 . ... 0.1.. PC . roptip. i, tb. protrcp. bo. r.tb.r co... ip.o .b. p.o.ro. 
vA* input data. 


Tb. 0, tb. POP.,.. mw (I, a, p, z, 22, 

tor X .... .. .. ipdopopd,., ^ 

tl... Of 1. no. . 2 ^ ^ 

.. tb. ..rrepopd,., X. a. .rr.p a 

IPdopopd... o.ri.b,.. X ,p . op.-di,.p.i..., ,rr.p. a. op.-dt....i...i .rr.p. 

n ppd ZZ . 0 ..,,. ,p , 4 . „b„ d.p,.d..t v.rl.bl... Si... 

...b Of a... .rr.p. ,a. XT. Z2) .op..,.. 

P.r«...r H i. ...d to t.di..t. bo. „pp ...... ... 

tbp. bo. lug, 0 . 0 b .rr.p 1. to b. dl....iop.d 1. TOLIN. 

I. obopld b. «,t.d tb.t tbi. roptlp. ... ..,2 b. ...d ,b.p tb. ...... ,. 

.rr.p a pro Up.. ,. ..pd..r...ip, a. ....!,. .f ,b,. a...r i.t.rp.- 

I...O. ...tip. .r. . 0 . ,.ud if tt. .Pi... a ^ 
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•r« Mt llsta4 in inoratsiat ordar* An aaanpla at tha aaa of tktia roatiaa oe~* 
enra vban a Talna of altitoda (X) ia input* Tha rontina will »at«ra valnaa of 
aabiant taaparatua (f) and praaanra (Z) Tia llaaar intarpolation fxoai tha ta»- 
paratnra and praaista arraya* TY and ZZ. 

Upon antarini tha rontina* a aback ia made to datarmina if tba indapandant 
variabla I ia laaa than tba firat (aaallaat) valna ia tba array II. If it ia* 
than tba firat clanaat of array YY and tba firat alcBcnt of array ZZ ara ra- 
tnrnad aa T and Z« raapactivaly* to tba ealliag rontina. If tba indapandant 
▼ariabla X it graater than tba laat (largaat) alaaent in array XX* tban tba laat 
alaaant of array YY and tba laat alaaiant of array ZZ ara ratnrnad aa Y and Z* 
raapaotivaly* to tba calling rontina. If tba valua of tba indapandant variabla 
X ia not out of tba range of tba array XX* tban linear intarpolation ia aaployad 
to datarsiaa tba aorraaponding valnaa of Y and Z* 

2.5 8UBROOT1NB TINTd 

Snbrontina TINTS ia a linear interpolation rontina. Tbia rontina intarpo~ 
tlx ot laaa dapandant variablaa aianltanaonaly aorraaponding to an input 
indapandant variabla. Tba argnaant of TINTS eontaina IS alaaanta, Tba firat 
two alaaanta ia tba argnaant ara valnaa of tba indapandant variabla and tba 
array of atorad valnaa of tba indapandant variabla. Tba neat 12 loeationa of 
tba anbrontina argnaant oeanr ia aim paira aorraaponding to tba aix dapandant 
variablaa. Tba firat alaaant of tba pair ia tba value of tba dapandaat variabla 
ratnrnad by tba prograa* and tba aaeond alaaant of tba pair ia tba array of 
input valnaa aorraaponding to tba valnaa in tba indapandant variabla array. Tha 
indapandant and dapandant variabla arraya nay aaob oontaia a aaaiana of 10 
valnaa. Tha laat alaaant of tba anbrontina argnaant aorraaponda to tba anabrr 
of valnaa aontainad in tha indapandant variable input array (10 or laaa). 


Tb« alaMaot of tba •nbrontina argnnant eoataiaa as intagar 
>ataa. It tha valaa of tbia Intagaa is graatar than 0, than TINTtf parforat In- 
tarpolatloa aatwiag tha iadapaadaat varlabla la aoaiaoratalag. If tba aalaa of 
tbla latagar it lata tbaa taro« tbaa TXNT6 parforaa llaaar iatarpolatioa tttaai'- 
iag tba iadapaadaat variabla to ba aoadaeraatlag. 


Now eoaaidar tba eata wbara tba iadapaadaat variabla it ooadacraaaiag. 
Bafora iatarpolatioa it attaaptad* tba iapat walna of tba iadapaadaat variabla 
it ebaekad to dataralaa if It la latt tbaa tba first aleaeat la tba iadapaadaat 
variabla array. If tba input value it latt tbaa tbit first alaaaat, tbaa tba 
kba dapaadaat variablat raturnad to tba oalling routiaa ara tat oqual 
to tba first alaaaat of ateb array. If tba iadapaadaat variabla it largar tbaa 
tba first alaaaat* tbaa aaotbar ebaok occur t to dataraiaa if tba iadapaadaot 
variabla it graatar tbaa tba latt alaaaat ia tba iadapaadaat variabla array. If 
it it* tbaa tba valuat of tba dapaadaat variablat raturaad to tba calliag 
routiaa ara tba latt valuat ia aacb dapaadaat variabla array* If tba iadapaa~ 
daat variabla pattat both tattt* tbaa TINT6 intpactt to ratolva wbat^.ar t'x aot 
tba iadapaadaat variabla it aqutl to oat of tba valuat ia tba iadapaadaat vari- 
able input array. If equality occurs* tbaa tba value returned fox aacb dapan- 
dant variabla it tbat value ia aacb dapandant array corratpoadiag to tba inda- 
paadant variabla. 

If tba value of tbo indapendant variabla it not found to ba latt than tba 
^^*tt alaaaat* graatar tbaa tba latt alaaaat* or aqutl to any of tba alaaeatt ia 
tba input array* tbaa linear iatarpolatioa occurs using tba two alaaanta of aacb 
dapaadaat array wbicb corratpondt to tba two input iadapaadaat variablat brack- 
eting tba input iadapaadaat variable* 


An nnalofoot proc«4nr« it ptrlocntd by T1NT6 if tbt indtptndtttt varitbln it 


noninorttting intt«itd of &oadt«r«ttini. 


Tbit roBtlnt it written for oitbtr t nondtvrtttini indtptndtnt wtritblt or 
• aoainertttini indtptndtnt wtritblt* toch tt tint inertttiag dnring t trtJtoto~ 
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2,6 FUNCTIONS INTPl AND INTP2 

Tbttt fnnttiont art nttd to ptrfora liattr inttrpolntiont in ttni 

lof^^ tptet. Tht nrtnatntt of INTPl trt: 

I ” (input? Prtttnrt iadtn 
J2 ” (inpnt) Bifbtat ttaptmtort index 
J1 - (input) Lowett ttnperttnro index 
K - (input) Tbemtl property index 

P “ (input) Deltt wtlne of independent wtritble it (X~X1)/(X2-X1) 

Vtlntt tre interpolated fron tbe 6(1* 7* E) nttrix* 


Tbe arguientt of INTP2 tre: 

X “ (inpnt) Lowett wtlne of dependent wtritble 
T - (inpnt) Bigbett wtlne of dependent wtritble 
F - (input) Deltt wtlne of independent wtritble 

Both fnnetiont linearly interpolate in loglO tptoe of tbe dependent wtritble* 


( START ^ OWCilNAL PACV ; .. 

Vs^ J OF POOR QUALITY 

r ZERO ARRAYS I 

r 

I READ INPUT DATA. W(H ARrFI 
f WRITE OUT INPUT DATA 1 

r 

f CHECK INPUT DATA FOR ERRORS I 

T 

I'lNITIALIZE PARAMETERS I 

r 

I BE.IN CALCULATION LOOP I 

— - ■! 

I COMPUTE FREE STREAM PROPERTIES 1 

Lcompute boundary layer edge conditions I 

i — * 

( COMPUTE CROSS FLOW CORRECTED RUNNING LENGTH I 
r VIRTUAL ORIGIN CORRECTION FOR TURBULENT FLOW | 

I — ^ 

1 COMPUTE LAMINAR AND TURBULENT HEAT TRANSFER COEFFICIENT ! 

A ' ' 

I RADIATION EQUILIBRIUM TEMPERATURE I 


r WRITE OUTPUT I 

I 1 



Fig. 2.1 Routine MAIN 




Ssotion 3.0 


FLUID PSOPEtn R0DT1NB8 

Ir tliis ••otlon. th* options of eospatini tho Irso'-ttros*.'' sir propsrtiss 
sro sot forth. Thoss iaolndo idssl sod rssl fss offsets. 

3.1 8UBR0DTINB AH62 

This is oas of th# otaosphsr# routines in NlNIVfiR. The rontiae deseribes 
the U. S. Stsndsrd Atnosphere. 1962 (ref. 2). Oivsn altitude (ft>. the rontiae 
yields eabieat teaporeture (degrees R). density (slng/ft*). pressure (Ibf/ft*). 
end speed of sound (ft/seo). The altitude rsage of the table is froa 0 to 
2*068.776 feet sboye see level. The theraodynsaie propei‘*’i«e. toa'>err.ture end 
density* ere exponential eurve-fits to the U. 8. Stsndsrd Ataosph. <(, 1962. The 
eabieat pressure end speed of sound are eoaputed froa vK<. i^' '> . equations. 

This routine does not esll say other rout.in> xns only to the ea^ling 

routine. 

9.2 8DBR0DT1NB AXMS4 

This is the 1963 Petrieh etaosphere routine. The input into the routine is 
altitude (ft) end the routine returns aabieat teaperature (degrees R)* density 
(slug/ft*)* pressure (Ibf/ft*)* and speed of sound (ft/see). The range of alti- 
tudes is froa 0 to 2*300*000 feet with errors of less than 3 pereent below 
400*000 feet. This routine is entire within itself and does not oall any other 
aubroutines. The routine breaks up the altitude span into aeveral rangea. Baeh 
altitude range then has the ataospherie properties ealeulated via a coabination 
of polynoaiala and exponential eurve-fits. 

9.3 8DBRODTINB VRA71 

This is the Vandenberg 1971 reference etaosphere. The input into the 
routine ie sltitude (ft) and the routine returns aabieni teaperature (degreea 
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R), . ity (tlat/ft*), prttsnc* (Ibf/ft*), and tp««d of toand (ft/teo). Tbo 
rtofo of altitado it fr«B 0 to 2*300,000 ft. This roatiao does aot call aay 
otbar roatlaaa. Tba aabroatiaa ealoalataa aablaat viaeoaity aad aolooalar 
«al|ht, bat tbaaa aalata ara aot aaad. 

3.4 ROUTINES HOLIER. BLOCE. BLOCEA AND SLOPE 

Hioao ara tb# roatlaea ooaatitatiat tba air Molllar diagraa and act aa an 
aqaatloa of atata to accoaat for raal gaa affaota. Tba data ara eontalnad in 
tba 0(1. 7. E) aatriz. 


0 ( 1 . J. 1 ) 
0(1. J, 2) 
0(1. J. 3) 
0(1. J. 4) 
0(1, J, 5) 
0(1. J. d) 
0(1. J, 7) 


- Taaparatara 

- Enthalpy 

- Entropy 

~ Viacoaity 

- Prandtl Naabar 

- Spead of Soand 

- Danaity 


BLOCK 

BLOCEA 


Tba bloeb data BLOCK containa tba praaaara array PT(14) and tba variablaa 0(1. 


J, K), K - 1,4. Tba block data BLOCKA containa tba variablaa 0(1, J, E), E » 


5,7. 


Foar optiona may ba aaad to lo.^k ap propartiaa baaad on two independent 
variablaa. Tbaaa optiona ara: 

Non “0 P.H (Praaaara. Enthalpy) 

1 P.S (Praaaara. Entropy) 

2 P.T (Praaaara. Tamparatara) 

3 B.S (Enthalpy, Entropy) 

Baaad on the optiona aalactad. tba limita of tba program ara ebaokad 

10“‘* i P i 10* Atm. 

48. 1 B i 75580 Bta/lbm 
200. 1 T i 30000 R 
.525 i S/R 1 12.368 

If tba indapandaat variablaa ara oataida of tba preceding range of validity, the 
roatina aata all tbarmodynamic variablaa to aero and rataraa with HHBAL • 1. 


Oivan tba indapandant variablaa. tba dependant variablaa ara obtained by 
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interpolation in tbo 6 natris. For NOFT " 0, 1 , and 2« tba intarpolation pro<> 
cadnra la tba aaaa. Tba praaanra indas la found corraapondlni to tba table 
praaanraa below and above tba Input praasura level* At tba lover praaaura 
level* tba dependent varlablea are detemlned by two ealla* Plrat* tba funetlon 
SLOPE la ealled to oaleulate 

I - XI 

Fn • 

12 - XI 

where X la tba loglO of tbe other independent variable and the XI and X2 value a 
are tbe loglO valuer of tbe neareat pointa in tba 6 aatrix for tbe lover prea** 
aura* Given tbe value of FTl* all of tbe dependent varlablea are oaleulated by 
linear interpolation in loglO apace uaing FTl in a call to INTPl* Libewiae* tbe 
dependent varlablea are detemlned at tbe higher preaaure level in tbe tablea* 
Finally* tbe dependent varlablea are detemlned for tbe input preaaure level by 
linear interpolation in loglO (preaaure) apace by calling INTP2* 

For NOFT “ 3 tbe interpolation procedure ia aiailar but uaea different 
logic* Tbe 0 autrix ia aorted to obtain tbe aurrounding pointa of enthalpy and 
entropy* Tbe nondimenaional entropy ia uaed in a call to SLOPE to obtain tbe 
paraseter 



at tbe lower preeaure level where 8 ie tbe loglO of tbe entropy* A call ia made 
to INTPl to obtain tbe dependent varlablea at tbia preaaure level* Tbie proceaa 
la repeated for tbe higher preaaure level* Finally* tbe dependent varlablea are 
detemlned at tba given enthalpy level by a call to INTP2* 


Aft*r the oalcnlttlos* «&d iaterpolatioai azf asdt, soae 7 «rlabl«s are coe- 

vextad la naits aad tha coapratalbility aad apaelfio hast ratio ara aalcelatad 
for output . 

Tha arguaaata of tba call for HOLIER ara 

H - (laput/output) Eatbalpy 

P ~ (laput/output^ Praaaura 

NOPT - (iaput) ladapaadaat Variabla Optloa 

T - (laput/output) Taaparatura 

21 - (output) Coapraa Sibil Itp Factor 

S “ (output) Batropp 

R - (output) Daaslty 

G1 - (output) Spaeifle Haat Ratio 

la additioa tba ▼iseosity la output via tba eoaaoa statanaat /BANE/. 

3.S SDBROBTINB BINIRP 

Subroutiaa BINIRP eaa ba uaad to ealeulata tba Prandtl auabar Pr, Lavis 
auatbar La* and diaaooiatioa aatbalpy Ed for air. Tba Prandtl aumbar is coaputad 
as a fuaetioa of tanparatura (dagraas R) and as a function of tba ratio of tba 
prassura at tba atagnatioa point to standard praaaura. Tba Lavis nuaibar aad 
diaaooiatioa aatbalpy ara oaloulatad as a function of tanparatura (darraaa E) 
aad as a fuaetioa of tba ratio of tba dansity at tba stagnatioo poiat to tba 
dansity at atandard conditions (2,37 * 10 * slug/ft*). Routine BINTRP is oon- 
poaad of tbraa tables. Table I corraspoada to dissoeiatioa aatbalpy* Table II 
aorrasponds to Louis aumbar* aad Table III corresponds to Prandtl aunbar. Tba 
first too alamaata of tba routine argument ara tba tanparatura (dagraaa B) aad 
tba * assure ratio* or tba tanparatura (dagraaa R) aad tba daaslty ratio. Tba 
third alanaat of tba routlaa argunant is tba table aunbar (I* II* III), Tba 
last (fourth) alanaat of t* routine argunant is aitbar Bd* La* or Pr, dapaadiag 
upon tba input. Upon antaring this routlaa* various ebacks occur to insure that 
tba input values ara not lass than those values for which a calculation can ba 
parfornad. If tba tanparatura is nagativa* tha pressure ratio la lass than 10*”* 
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or tho 4«ntlty ratio it Ittt than 10 * tbaa an arror naaaafa ia priatad and 

NEBROR ■ 1 ia rataraad to tba oailiag roatina* If tka taata art paaaad« than 
BWIRP proaaada to aaa doabla iatarpolatioa to dataraiaa tka appropriata oatpat. 
If Bd < 0, tkaa Bd • 0 ia rataraad. If Bd > 12450 Bta/lba. tkaa H4 - 12450 
Bta/lba ia rataraad. If La < 0.5, tkaa La - 0.5 ia rataraad to tka aalliat 
roatiaa. If Pr < 0.678* tkaa Pr • 0.678 ia rataraad. If Pr > 0.754 and 
T > 24000 datraaa B« tkaa Pr ■ 0.754 ia rataraad. Praaaatly tkia roatiaa ia 

aaad only to datarsiaa Bd and Pr. Tkaa, tka otkar oaa>tkird of tkia roatiaa ia 
aot aaad. 

5.6 SDKtOTTlNB HANSEN 

Sabroatiaa HANSEN ia eallad by roatiaa FLOW to ealealata aiaeoaity aa a 
faaotion of taaparatora and prataara. Tkia aabroatiaa oallt HOLIER and rataraa 
vlaoo.ity Tia ooMoa HANK. If tka yi.ooaity i. pr.atar tkaa tka awitek yalaa of 
0.29551E-08 alag/ft aao* tka viaooaity from HOLIER ia rataraad. If tka yiaooai- 
ty fro. HOLIER ia laaa tkaa tka awitok valaa, tka yiacoaity ia oaloalatad baaad 
on tka ralationa of Haaaaa and Halna (NACA IN 4559). 

9 ■ B/D 

vkara 

IL - lof^^ P/2116. 

A - ((T/1800) (1-0.125IU - 6.5)/(1.5 ♦ 0.125R) 

B • 1. ♦ 0.025 (T/1800. )<1 ♦ Taak(A)) 

C - (T/1800 - 14.5 - 1.5a)/ (0.9 + o.ia) 

B - a® ♦ 1.0 

9, - 2.27 a lO"* T».*/(T ♦ 198.6) (Satkarlaada) 

T ■ Taaparatara (Dagraaa R) 

B “ Praaaa^'a (Ibf/ft*) 

If tka yalaa of B/D < 0.04 tkaa B/D ia aat aqaal to 0.04. 
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Section 4.0 
FLOWFIELD ROUTINBS 


me teetlon of tbo report deals with those snhrontines that are need to 
ooapnte the local flowfield and boundary- layer edge conditions around the body 
under consideration. The nethods used to evaluate shock angles* flowfielda 
•hocks, and Prandtl-Meyer espansions are presented. 

4.1 SOBROOTINE FLOW 

Subroutine FLOW is called by MAIN. Ihe subroutine argunent contains two 
paraaeters* FF and ALPHA. The element FF is a one-dimensional array containing 
a possibility of nine consecutive flowfield and pressure flags (always occurring 
in pairs with the exception of the swept cylinder flowfield FF « 39 and the 
Prandtl-Meyer expansion flag FF - 29). The element ALPHA is a one-dimensional 
array containing a possibility of nine consecutive flowfield angles correspond- 
ing to the flowfield flags. The transfer of other information, such aa upstream 
conditions, downstream conditions, etc., is performed through common blocks. 
The free-stream temperature and pressure are first needed to compote f^«e-stream 
viscosity via subroutine HANSEN. Ihe first check performed by FLOW oeterminea 
If the Mach number upstream M^ of the first flowfield flag la greater than 
«alty. If M^ < 1 and thie is the firsw or second flag in the array FF, then the 
program ignores the flag and proceeds to calculate stagnation conditions baaed 
on the isentropic ideal gas relations. If the flag following the M <1 verifi- 
cation calla for a Prandtl-Meyer expansion (FF - 29), then a deceleration of the 

flow occure. The flow decelerates to a higher pressure calculated by modified 
Newtonian theory. 

tnowing the local stagnation conditions and a local subsonic modified 


14 


N«wtoal«a prattos** *11 of tho otbor local pxopaatlaa ax# dataxalnad by tho 
ideal taa iaoAtxopio xolatloaa. If tha daoolaxatioo a&|lo la gxoatax tbaa ox 
# 4 «al to 90 datxaoa* than tha local aondltioaa axe act agnal to ataiaation aon- 
ditioaa* Pox anbaoala flow* tha local aonditiona dataxniaad abova axe alao act 
aqnal to tho local aonditiona along tho atagnation lino of an infinito anopt 
aylinda«. laving tho local flovfiold aonditiona, tho Eehoxt xofoxonao pxopox- 
tioa axo nozt aalenlatod. If tho hoot txanafox flag ia ognal to 5 ox 7. than 
tho PjpHjp pxopoxtioa axo alao aalenlatod, Tho xontiao than xotnxaa to MAIN, 

If M <1 and thla ia not tho fixat flowfiold ox pxooanxo flag, than tho 
H 

local atagnation conditions axo takon to bo tho saao as thoao aoxxoapondiag to 
tho aonditiona downatxoam of tho pxoviona pxooanxo ox Pxandtl-Moyoi flag. If 
tho flovfiold flag following 1^ < 1 «oHo fox a Pxandtl-Moyox ozpansion (FF - 
29), than a doeoloxation of tho flow oecnxs. Tho flow docoloxatos to a highox 
pxosanxo again aompntod via nodifiod Newtonian thooxy. Knowing tho pxossnxo and 
also hawing tho ontxopy fxon tho (ieontxopie) npatxoan conditions, tho othox 
local thoxBOdynania pxopoxtioa axo dotexaiinod by nsing tho Molliox diagxaa (an— 
bxontino HOLIEK). If HOLIER xotnxns IDEAL > 0, than tho pxogxam svitahoa to tho 
idoal gas aalanlation pxowionaly diaansaod fox < 1,0, If HOLIER xotnxaa 
IDEAL • 0, than FLOW ahoaka to dotoxnino if tho doeoloxation aaglo ia gxoatax 
than ox oqnal to 90 dogxooa. If tho doeoloxation aaglo is gxoatax than 90 do- 
gxoos, than tho local conditions axo sot oqnal to tho local atagnation aoadi- 
tioaa. Fox snbaoaie flow, tho final eonpntod local conditions axo alao sot 
oqnal to local conditions along tho atagnation lino of an infinito awopt cyl- 
indox. Knowing tho local flow pxopoxtioa, tho Eekoxt xofoxoaeo pxopoxtioa axo 
aost aalenlatod. If tho host txanafox flag is oqnal to oithox 5 ox 7, than tho 
9,9, pxopoxtioa axo alao aalenlatod. Tho FLOW xotnxaa to HAIN. Tho diffoxonco 
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b«tw««n < 1 aod FF - J or 2. and < 1 and FF > 2 la tbat for FF - 1 or 2 
tha flnld la traatad aa an Idaal gaa, and for FF > 2 tba fluid la traated aa a 
raal gaa. Only tba aqnatlona of atata ara dlffarant. 

Now If > 1, FLOW ahaaka to dataralna If FF * 35» FF • 36, FF ■ 38 (thara 
la no FF • 37) or FF ■ 39. Each of tkaaa flowflald flaga corraaponda to a dlf- 
farant t< cinl^ua for dataralnlng tka akoak angla fi. For FF ■ 35 or 36, tka 
▼alua of ALFA la takan aa tka flow daflaatlon angla for a akarp wadga or akarp 
aona, raapaatlwaly. For FF - 38, tka akock angla la raad Into tka progr.« aa 
ALFBA. .And FF ■ 39 oorraaponda to tka awapt eyllndar flowflald and tka akoek 
angla la aqnal to ALFA. Sukroutlna FLOW (for FF • 35 or 36) calla routlna PCS! 
to dataralna tka attackad akoek angla aa a fnnetlon of H and tka wadga or eona 
daflaotlon angla. If NESROR 0 la raturnad from PCSW, tkan tka input daflac- 
tlon angla auat kava baan altkar nagativa or aqual to saro. FLOW tkan Immadl- 
ataly raturna to MAIN and tranafara to tka nazt data eaaa. Tka naxt akaek la to 
datarmlna tka walua of OFT. If OFT - 1, tkan tka flow upatraam of tka wadga 
muat kava baan aubaonle. If OFT ■ 2, tkan tka akock aaaoclatad wltb tka wadga 
la dataokad and tka akoek angla la, tkarafora, takan to ba 90 dogr'saa. (If tka 
akoek la dataekad, FLOW parforma a normal akoek azpanalon (modified Newtonian) 
to tka praaaura flag angla aaaoelatad wltk tka flowflald angla.) If OFT - 0, 
tkan PCSW kaa aueeaaafully datarmlnad tka akoek angla aa a function of upstream 
Mack number and daflaotlon angla. If OFT - 1, tkan tka routine proeaada to make 
a aubaonle ealeulation aa dlaeueaad above. 

Hawing obtained a akoek angla fi (altkar Input for FF • 38, p • ALFA for FF 
■ 39, or for a akarp wadga or eona from PCSW) (90 dagraaa for a dataekad wadga 
or eona akoek), FLOW next ealla DWNSTM to datarmlna tka propartlaa downatraam of 


th« thoek. Sokrotttln* OWNSTN «••• p aloof with aw air Molllar 4iatraa (a»- 
baootlaa NOLIBR) aa an aquatioa of atata and tha apstaaaa aoaditioaa to aoapota 
tha eonditioaa dowaatraaa of tha ahook* If tha propartlaa raqolrad fcoa tha 
diagraa aaaaot ha ohtaiaad* than IDEAL ■ 1 or 2 Bay ha ratoraad to FLOW* 
If HEAL ■ 1« than tha rootlaa proeaada hy aaplbyisg Idaal gaa ralatloaahlpa 
(dlaeaaaad latar). If IDEAL • 0 or 2* than tha prograa proeaada hy aaployiag 
tha Kolliar diagraai aa tha aqwatloa of atata. 

Oaea tha eoaditiona dowaatreaa of tha ahoek are kaowa aad aaaoalag IDEAL * 
0. FLOW thaa proeaada to datanaiaa tha local wall praaaura haaad oa tha praatara 
flag (14. 15, 16, 17, 18) followiag tha flowfiald flag (35, 36, 38). If PF ■ 
Sahroatiaa TBLIN ia eallad to iatarpolata tha praaaura ooaffieiaat froa 
aa laput aat of tabular waluaa with ooaffieiaat of praaaura giwaa aa a fuaetioa 
of fraa**atraaa Naeh auaihar. If FF • 15 or 16, thaa auhroutiaa PCS! ia eallad hy 
FLOW to datamiaa tha praaaura ooaffioiaat for a aharp wedge or ooaa aa a fuao- 
tloa of aad daflootioa aagla. Agaia if NERSOI ^ 0 ia roturaad froa PC8W, 
thaa FLOW raturaa to MAIN aad akipa to tha aazt data eaaa. If OFT • 1 la ra- 
turaed, thaa tha prograai awitehaa to parfora a auhaoaio oaloulatioa. If OFT ■ 2 
ia raturaad (dataohad ahook), thaa FLOW eoaputaa tha local aodifiad Nawtoaiaa 
praaaura. If OFT • 0 ia raturaad, thaa tha loeal edge of tha houadary-layar 
praaaura ia dataraiaad haaad oa tha praaaura ooaffieiaat raturaad froa PCSW. If 
FF - 18, thaa tha local praaaura ia dataraiaad froa aodifiad Nawtoaiaa theory. 
Tha other loeal tharaodyaaaio propartiaa arc ohtaiaad froa tha Mol liar diagraa 
routiaa (NOLIER) aa a fuaetioa of local edge praaaura aad tha dowaatraaa aatro* 
py. taowiag tha loeal aathalpy aad tha atagaatioa aathalpy, tha walooity ia 
coaputad froa tha aaargy aquatioa. 
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urr. 17. tk. 4. a,. 

.,,1,. n„ ^ 

.»... rt..k ..„. 1. „ „„ 

•tlo# •■|U. Ik* *keek tail* ■*•* b* knows *ad ispwt nok*r PP - 3$, ji* 

.«|. *f tk. b...d.w3-U3„ „l„,„ 4 . „„„,4 4 „, ,„f.„ 

»|l«-.f-.tu«k. tk. rt.,k ..,1. a, .,.tr... T.l..itp. Ik.. ,k. 1...I ... 
tk.l„ 1. *«,.t.4 l.„ u.. .„w,3 .,„44.., c.,4., ,4. 

k.ki.4 tk. rt..k tr« OINSn (PP . 33), tk. .tk.t I...1 tk.,...r„. 4 . 

pt.p«tl.. „. ,bt.4»4 .4. a, ,i. «.m.w 41.,.„ „.,4„. „„ „,j 

U .. ..tp.t ot tk. 11.114., dl.„„, 43., 3.4., ,34.4„^ ^ 

.». «.fft.l..t („» PCS,, .. ... 3.„ „ . „ „ „ 

..3l«f-.tt..k 1 . ,.„., 43« tk. rt..k ... 1 ., . 3 .. 43 . 4 . 

««P.t. tk. .4,. «,.4ttl... .4. tk. p.,.11.1 .30.3 t..k.t... (,4tk p .. tk. 

.k.,k ..,1.,. 14 4kU 43.. ..y 1.11„4., ll„fi,43 44.,. .„ 4„„,.^ 

•ad the f€ogT9m proceeds to eonpato o it md/ar » 

»p»*c saa/or Eckert refereace properties be- 
fore retaralag to MAIN. 


M PP - 14, IS. 14, 17, IS ..4 .. .k.o„.144y ...... ..,3 ,. , > 3 (p, . 

14, . tk.. k.4„. p,....44., 4. „„... 43. „„ ^ 

«.t.»l» 44 . P,..4tl-lb3., ..P...1.. 1. 4..U.4, 14 . P,..4tHU7., ......1.. 

t. 4..1..4, tk.. PM, ..11. BMW 4. ..,..,.4. ...4444„. 4.,..4„., 43, 

pt..).., 14 « P,..4tl-lb7., ..p..tl.. 1. 4.tl,.4, tk.. nm bpp..,.. PKEIW. 

H..t . .k..k ...„. 4. 4.,.„1« 44 .441.1..., ,.„41.,4 4,.,. (,!» .. ,4.3... 

... „.„.4. 14 43.,, ,3„ 43, ...41.M. 1..., 

M.41tl.u .,* ttkt. .. sp.t,.u oo.41tl.a. 4 m tk. uit 41w41.14, 14 .0 

nMll.14 n.„ ... ,,....4, 43„ P3„ ^44,444^4 44 ..1..1.4, ,4.4 b.3.,4 

tsfsreaee properties. 
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f-'. 

W FF • 39 it nttd, tbtn FLOW sttt tlit downtirtts t&tsopy tad tbt tttgat- 
to «o«putt tbt othtr local tbaraodyatale properties. HsTlaf 
tbe local tberaodyaMic properties, tbe local ataiaatloo sad static flow proper* 
ties are coapated. Tbea tbe ead/or Eckert refereace properties are eoapat* 

ed. Ibe swept cyllader flowfleld optloa asp be used aloae or as a laat flow- 

field flag ia a series. Aay flc .'laid optloa atsd after FF - 39 will be Igaored 
by FLOW. 

If ia tbe ase of tbe Molller diagraa roatiae tbe error aessage IDEAL • 1 or 
3 is retaraed to FLOW, tbea FLOW skips tbe real gas ealcalatioa sad perforas aa 
••• oalcalatloa. Ibe Nollier diagraa roatiae (HOLIER) is capable of baa- 
dllag botb ideal aad real gas calcalatioas. The ideal gas relatloas ased la tbe 
prograa are priaarily ased wbea aotblag else works. Ibis is to keep tbe prograa 
;i| raaaiag rather tbaa stop the prograa becaase tbe Nollier diagraa roatiae eaaaot 

extrapolate to a valae off tbe table. Oace FLOW awitebea to tbe ideal gas cal- 
j calatloa, all tbe flowfield aad pressare optioas are exeeated asiag tbe ideal 

j gas relatioasbips. Tbe ideal gas eqaatloa of state is iaplicitly ased iastead 

of HOLIER. Tbe shock eagles aad pressare coefflcieats are deteraiaed exactly 
,4 the sane way as for a real gas ealcalatioa. Ibe dowastreaa eoaditioaa are aow 

^ deteraiaed by tbe ideal gas roatiae DOWNID iastead of DWNSIN aad tbe eoaditioaa 

^®*®*****“ P»*a4tl-Meyer expaasioa are coapated ia tbe ideal gas roatiae 

^ PMIJ) iastead of PNBXFN. All logic aad error checks are precisely tbe saae as 

for tbe real gaa ealcalatioa iaeladiag tbe ideal gas deteraiaatioa of tbe o a 
ead/or Eckert refereaee properties. 

4.2 SUBROOTINE DWNSIN 

Sabroatiae DWNSIN is eslled by FLOW to ealealate tbe eoaditioas dowastreaa 
•» . .koekw... .. 
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th« downttr*** oonditioBt tot both r«al and idaal laaaa. Bowtvar, tha calcala" 
tion proefdnra la baaad upon tha aatomplioa that a raal gaa eaicalatioB vill 
taka plaea. If aabroatloa HOLIER daterniaea that aa idaal gaa ealealatioa la 
aaffieiaat* thaa tha idaal gaa aquatioa of atata ia aaad iaataad of tha Mol liar 

diagraa Talaaa* 

Tha ahoak aagla fi ia iapat iato DVNSTM through tha aabroutiaa arguaaat. 
For aa attaehad vadga or eoaa abock, R ia dataraiaad ia aabroutiaa PCSW, If tha 
ahook ia dataraiaad to ba dataohad by PCSW* thaa fi it takaa aa 90 dagraaa for 
iaput iato DWNSTM. If tha parallal shook optics it ohotaa* thaa tha shock aagla 
ia aat aqual to tha affaetiva aagla-of-attaek. 

Haviag aa iaput shook aagla p DWNSTM aazt requiras tha iaput of eoaditioas 
upetraaa of tha shook. Tha upstraaa eoaditioas are iaput iato DWNSTM via a cob- 
Boa bloek with FLOW. Oaoa tha upstroaB eoaditioas sad shookwava aagla ara 
*^*^^*^1*» DWNSTM prooaads to ealeulata tha dowastraaB eoaditioas. First* tha 
aoraal ooapoaaata of tha upstrasB valoolty* aad Msoh auaber ara ociouTatad. 
Thaa aa iaitial guess It aada for p,/p^. If 1 4000 ft/sac* tha ideal gas 
aquatioaa ara used to aaka tha iaitial guess. If > 4000 ft/aac* tha iaitial 
guaaa for p^ la aada usiag various aapirloal aquatioas dapandiag oa vhathar 4000 
ft/aac < Vj i 8000 ft/sae* 8000 ft/aeo < i 14000 ft/aao. J4000 
ft/aao < Vj i 26000 ft/sae* or 2600U ft/sae < V^. Thaa a ohaok oeoura to iasura 
that Pg/pi > 1*0* Subroutiaa DWNSTM thaa bagiaa tha first of a Baziaua of SO 
tt 0 tttlou$ to dataraiaa tha dovastraaa propartiaa. 

a 

This diseussioa aaploja tha followiag aubseripts: 1* 2* N* T. Tha subserlpts 
1 aad 2 oorraspoad to eoaditioas upstraaa aad dowastraaa of tha shook* raspae- 
tivaly. Subserlpts N aad T eorraspoad to velocity eoapoaeata aorBsl aad taagaat 
to tha shook* raapaotivaly. 
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Bavint an Initial gnaaa for tfaa 4ownatraan danaity tha downatcaaa 

nornal coaponant of yalooity, V,„, ia ealonlatad froa Equation (4-1). 




and than tha dounatraaa anthalpy and praaanra ara coapntad haaad on p, 

h. . h. * !!L 
2 2 
» m 

Pa ■ Pi ♦ [piVij,* - 

Now p.^ ^ ia datazainad froa anhrontina HOLIER via Equation (4-4). 

(a) 

Pa' - pChj. Pj) 


(4-2) 


(4-3) 


(4-4) 


« <*) (») 

„ < 0.01 

convarganca ia eonaidarad to haya baan 
aohlayad. If aonyarganaa i. not achiayad. than p.<‘> i. r.pln«,4 by and 

tha aalculationa. baginning with Equation (4-1) ara rapaatad. Thia prooaaa i. 
continuad until two oonaaoutiya itarationa ara within 1 parent of om anothar. 
®r until a naninun of 25 itarationa ia attainad. Aftar 25 itarationa altarnata 
danping eoaffioianta ara choaaii to obtain a guaaaad danaity 

- CP<‘> * (1 - o)p<‘> 

whara a - 0.7 batwaan 26 and 35 itarationa and • - 0.4 batwaan 36 and 50 it.ra- 
tiona. Onea aonyarganaa ia obtained, all of tha atatia oonditiona behind tha 
ahoah ara hnown yia Equationa (4-1). (4-2). (4-3) and anhrontina MOLIEB. Tha 
yaloaity behind tha ahoah ia given by 


Vi,^* 4 (Vicos3)‘ 


(4-5) 


21 


cot fi m y 


•T* 


*»u, ,k. .ondUl... b.M.d ,h. 
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j (4-6) 

DWNSl, d... „t ,.U.U.. u. ..k„ .......... .... 

«h. ,...11.1 .k..k (...„ 

P«.l.. ... ..l..l...d .f.„ ...„..., s.l..«... H.OT .... s ... j 

.. l.,«. .. »LIEH 1. d...„l..., .k. 

4.3 SUBROUTINE PCSW 

s.k....l« ,csw 1. k, .o„l„ no, .., ,.. 

t». .k..k ...I.. ..a ,....„. ^ ^ ^ 

roatia. cont.lat four table.: one for tbe ^ 

no for the ah.rp-wedge .hook .agio ,, , fwotloa 

OX upstxoftii Mach &iu&bcr and daf « 

..rt d.tl...l.. ..„., .k..p-..« .k,.V, ..„. 

.. . d....i« or 

.u.,-..d.. ,....„. .. . ^ ^ ^ ^ 

xlootloa .agio, .ad the l..t oao for th* 

the .h.rp-eoae pro. ..re ooeffieloat .. . 

foactioa of ap.tr.aa N.eh aoaber .ad doflootloa .a.l. r.o 

»oB .agio. Xao .rgaaeat of PCSW 

n„ ^ ^ ^ ^ ^ 

..d .k. .„o«i ,..„..., 

,..«>..., nABLB, 1 . .. ....,„ 

IldBLE - 1 UdU.,.. .k,.k ,.„, t„ . .k.r, ,.d,. 1. d..l„d 
ITttU - k Udlc... .k.,k ...1. f.. . ., 

»«U . S IkdU.,.. 

IMU - 4 Indlc... „....„ t„ . .k.„ ,.., 1. a...„a 
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Th« foortb OFFJBL* if flflcDiftfd in PC8V and aiajr bava on# of tbraa 
poaaibla Talnaa, If OPPIBL - 1.0 ia vatiurnad. than tba npatraaa Maah ntabar and 
tba daflaation anpla oorraapond to conditiona prodncint a dataabad abock. If 
OFPIBL « 0.0« than PCSW baa anooaaafnlly eoaplatad ita fnnetioa. Tba laat 
pataaatac in tba avbrontina axi' aant ia tba aaloa of tba abeck angl* o* praaaora 
aoaffiaiaat for aitbar tba vadga or eoaa« dapandiag apoa wbieb ia raqaaatad froa 
tba aabxoatiaa. 

Baaiat aaeoaaafallp entered tbe routine, aaaeral ebaeka oeenr. If tba 
apatraaa Maeb naabar ia graatar tbaa 26. tba (apatreaa Maab aaabar) ia aat 
a«aal to 26. If < 1. than tba roatiaa xaiarna OFFTBL - l.If 6 (daflaation 
angle) ^ 0. tbaa PCSW pxinta oat an error aeaaaga and ratarna NESBfitt “ 1. Tbia 
eaaaaa PCSW to ratarn to tba MAIN prograa. Tba MAIN prograa tbaa will akip to 
read inpat data for tba nest oaaa if W(642) - 0.0. or will print tba atorad oat- 
pat ia roatiaa VANOUT and tbaa read input data for tba neat eaaa if W(642) 
f 0.0. If 6 > 60 dagraaa. tbon tba aboek ia taken to be dataebad for tba ooaa 
and OFFTBL ■> 2.0 ie ratarnad to tba ealling program. If 6 > 5S dagraaa. tbaa 
tba aboek ia eoaaidarad to ba dataebad for tba aadga and OFFTBL • 2.0 ia ra- 
taraad to tba aalliag roatiaa. 

;4axt. tba aalaa of ITABLB (1. 2. 3. 4) ia ebaekad to dataraiaa vbieb table 
ia to ba aaad. For aaeb aalaa of ITABLB aaaaral ebaeka oeear to dataraina if 
tba aboek ia dataebad. If tbaaa taata are paaaad. tbaa PCSW proeaade to ealea- 
lata aitbar a aboek eagle or a praaaara eoaffieiaat* 

Tba aadga end eoaa teblaa for aboek eagle end praaaara eoaffieiaat. ere 
rapreaaated bj earae-fita. Tba praaaara eoaffieiaat end aboek eagle era 
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csrvWitUd at a fnnotloa of Mach nimhcr aad daflaotloa aafla. Each of the 

four tablaa ooaaiata of acycral patchca* aaoh eoaacapoadiag to varioaa rcgioac 
for each table. 


After the yalne of ahoek aagle or preaaure eoefficicat la «»ilealated. a 
flaal cheek la Bade to detcralae If p 2 90 degreea or Cp > 1.8. Thla iapllea a 
detached chock aad OPFTBL ■ 2.0 ic retaraed to the ealllag prograa. 

4.4 SOBRODTXNE DOWNID 

Sabroutiae DOWNID is called by cabrontiae FLOW to ealealate the eoaditioas 
dovaatreaa of a chockwaye (aoraal. obll^ac. aad parallel). Thla roatiac oaly 
perforas ealealatioac applicable for aa ideal gas. The ahoek aagle 8 is iapat 
lato DOWNID throagh the sabroatiac argaacat. The apetreaa eoaditioas are lapat 
lato this roatiac yla a coaaoa block with FLOW, Siaee aa ideal gas calealatioa 
la belag aadertakea. there is ao aeed for aa iteratioa proeeta as ia DWNSTM. 
The calealatioa procedare begins with the calealatioa of the aoraal coaponent* 
of frce~streaa Mach aaaber 

M^n - Mj sin a 


f*‘®®***®“ •■ploy* the following cabaeripta: 1, 2. N. T. Ihe aabacripta 
1 end 2 correepoad to eoaditioaa apetreaa and downatreaa of the shock, reapeo- 
tiyely. Sabaeripts N aad T correspond to yclocity coapoaeata roraal aad taaaeat 
to the shock, rcapeetiyely. 


original page fa? 

OF POOR QUALITY 


Tht dow&ttrtsa pr«s«ar«, ttnptratnre tad Haoli aaabot givn by 


^ ^ 2YM|h* - (Y» 1) 

Pi (Y 4 1) 


(4-8) 


Ti P* (Y - 1) Mijj* ♦ 2 
Ti Pi (y ♦ 1) 



(4-9) 


(4-10] 


Tb« total prostuo bobiad tbo .bock, tbo .t.tic d.a.lty, local apccd of 


aoaad, dowaatrcaa velocity aad atatlc catbalpy aad tivea by 

, , r (y - 1) ,1 

Pm "Pi 1 ♦ — : M* 


•1 


Pa ■ 


RT* 


(4-n: 

(4-12; 


da • VyRTi 
Va • Mada 
hi ■ «24Ti 

aad tboa DOWNZD actaraa to FLOf* 
4.5 SUBROUTINB PNBXIN 


(4-13) 

(4-14) 

(4-15) 


Sabaoatiaa PNBXFN la tbo real gaa Paaadtl-Noyor capaaaloa aoatiao for air, 
Tbo lapat eoadltioaa to tba roatiac arc apatreaa atatto proportlaa aad atagaa- 

tlo. 1, ,1,, ......1,. ..,1. j. J 

«oav«,« M «rror to tbo ooot. 



ni« MUnutiti tfutlai 4««ct»ia| th> ml (ti tspistioa 1. 



( 4 - 16 ) 


witk P • ** k * n /J 

•ad a - a(k,«), t - « coast. 

wksre a - a(k.s) is a faaotloaal rspraseatatloa of tke Molllar ekart (Sabroutiaa 
HOLIER). 


Eaowiat th* valaa of • oaa desires to deteraiae eatkalpy. p, «t tke ead 
Of rke ..paasloa faa. Oaoe k is deteraiaed, all other properties are iaaediate- 
ly available. To deteraiae h. the above eqaatioa aost be solved for F. Ihi. is 
aekieved by eaployiag tke foartk-order Runge-Kutta aetkod. 

If 86 degrees < 0p„ < 103.2 degrees, tkea 6p„, is divided iato twelve e«aal 
.egaeats for aaaerloal iategratloa. If 68.8 degrees < 6p„ < 86 degrees, then 

V ••PMotod iato 10 e«aal dlvlsloas. If 51.6 < 6p^ < 68.C degrees, tkea 

V ioto eight equal latervals. If 34.4 degrees < e < 51.6 degrees, 

thea ep„, is separated iato sia equal segaeats. If 17.2 degree < < 34.4 de- 

8tooo. tkea Opfl, Is divided iato four equal latervals. Aad. if 0^ < 17.7 g,- 
****** tkea is divided iato two equal divisioas. Baviag sueeessfully deter- 
aiaed P at tke ead of tke eapaasioa faa, all of tke properties dowastreaa of tke 
fea are returaed to tke aaia prograa. 

4.6 80BR0DTINB PNID 

n>i. It tt. 14..1 Pr,o.tl-llt7t, »v..tu. „,u„. 
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anib«r. Tb« output oeosioto of loool Nteh nombor, toapototwo, psotturo, 4oaoi- 
tp* oatboipp# opoo4 of ooaa4 tad aoloolty aftoc tbo oxpaaoioa* Tho Haob aaaboc 
aftoc tho oxpaaaioa faa ia datarBlaod bp aaiag tbo Nawtoa^Sapbaoa itoratioa 
aatbod. A aaxisaa of 25 Itoratioaa ia alloaad. 


Seotioo 5 .C 


BBATINO RODTINBS 

». ..b».tiw. ...4 t. ..l..,.,. 
ieoB«trl«t tz* pze* 0 Bt*d la this section. 

5.1 SUBROUTINE FAIRID 

s.»r.„i„ wan. 1. 4.U.4 ..4 1. ...4 

«. „i„ .4 . a. 

BNa. a. ,.4„., ^ ^ ^ 

... 4 «.,.«, wa. i. 4. ..... 

a. FAaa u. b... ..4.44.4 ,. ...,„. 4... 

- «..44.4 4.„ b...... 4. .444.4.. .. b...4..4 4.4.. 4.« b...... b, .b. 

...boA .4 F.4 ..4 «44.U. a. b.«.4..4 l.F« «« «.b.4 b.. b... „4.44.4 .. 
..1...... ,.b...4. 41„ b...4« b4 . «>414i..t4.. .4 .b. „I..4.4 ...44..,. 


a. 41... ..b....4« ..U.4 »4 wniD 1. tEOItt. E...4 .. .b. ,.4.. .4 

a*. «b. ...... 4.« .E4IHB. . A...4..4 l-FT, ....41.4 .. 4... ..4...4.. 44.. 

ealcnletloa Is asde. 

JURE » 1 Bonaii^yY 

* b.»4..4 ay.. 44« ...,„.4.. P. 4 .. b...4b, ..4..4..4.n 4 . « 4 . b.M4 
ta ...b.4 .4 F.y ..4 ,44^44 .. a...ia4 1. ab4. 5.4. a. .,..,44. 41....4.. 
.i.b ...b.4,y 4. .b..4b.4 4.« BINISP .. . 4....4.. .4 ,b. 

P«rstaze beRlnd tbe nozasl sRoek and tke ratio of the stssnatlon denulty behind 

a. ,b..b .b. 4.b.4.y .. ....4..4 ,U..,b..4. w.44.4.... 14 MmOB 

, 0 4. ...„„4 by BINaF, a.. F*«a ...... 4. 

d... ..... fcT4b, .b..4»4 .b. b... «.„4.r ...444.4... ENCt, FATEID ...«,. 
.. U. ..444., „„4„. K, 4. ... 4....4,.4.. 



HIRE ■ 2 UMftUd Plfl, 

A flow •tatutioa point Roatiag ealeulatioa lo aodo by eolllag 
8IBBAT. per it Mt tqatl to 2,0 for ontpat dotorlptloa pvrpottt. Tbo boot 
trtatfor eooffleitat it output through tho Itbtltd eoaaoa /FRSMO/ by tht vtri*> 
tbit IT2L* 

If IRE » J Pf Molacultr Plow 

A fwt aoltoultr flow htttlug ctlculttioa it atdt for t plttt ptrptadioultr 
to tht flow dirtotioa by etlliag FHHEAT, PCT it tot tqutl to 10.0 for output 
doteriptioa purpottt. 


Out Itbtled COMMON htt beta added to FAIRID end it it otlltd /FREIIO/. Tht 
▼tritbltt trtatf erred through thit eteteaeat ere: 


IRE ~ (output) 


PCI 

STIL 

H 

Aim 

NONOQN 


AI2 

CRT 

C8I 

ZETA 

CHNI 

AZ2CIL 

XI 


(output) 

(output) 

(output) 

(output) 

(output) 


(output) 

(output) 


1 Boundary Layer 

2 Rarefied Flow 

3 Free Molecular Flow 

A wtriable tpeeifying the ealeulttion aethod uted 
Beet traaefer fila coefficient bated on enthalpy 
Beat tranafer fila coefficient baaed on teaperature 
Angle of attack plua local body angle 
Integer value of W(d4d) 

If NON CON > 0 a rarefied flow calculation can 
be aade if IRE > l 

If NON CON ^ 0 only boundary layer ealculatioaa 
can be aade (i«e« the prograa 
operate a identically to the 
unaodified MINIVER prograa) 
Rarefaction paraaeter ZeX* froa STHEAT 
Stanton nuaber froa STSBAT 


Not derived froa this aubroutine or aatoeiated callt. 


Sinee SWCIL otllt FA1RID« thia tubroutint thould bt ditontttd htrt* No 
aodifiettlont have been aade to SWCIL. Bowever. if W(646) - NONCON it input at 
greater than aero and SWCIL it ealltd. than a rtrefitd calculation could bt per- 
if IRE > 1. Thit aeihod it not rtcoaaended. 
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5.2 SUBROUTINE REGIME 


Thl« rotttino d*t«nila«t vbioli flow ro|lBO is opproprloto for bootl&g eoleo- 
lotions. Tbo Naeb noaber (AN), woloclty (UI), doatlty (RI), toaporataro (T), 
proooaro (P). aad ebaraotoristio loagtb (X) axo iapat tbroagb tba argaaaat. Tba 
aabroatiaa ealla DWNSTM to doteraiao tba poat-abock eoapraaaibillty factor. Tba 
flow ragiaa ia dataraiaad by tba ralatioaa glvaa ia Tabla 5.2 aad rataraa tba 
iatagar IRE. 

IRE “ 1 Boaadary layar (ooatiaaaa) 

“ 2 Rarafiad (traaaitioaal) 

“ 3 Fraa Molaealar 

Sabrontiaa REGIME la eallad by FAYRID. SPCSI aad SWC1L2. 

5.3 SUBROUTINE STHEAT 

Tbia aabroatiaa calcalataa tba baat traaafar eoaffioiaat* HTIL. to tba 
atagaatioa poiat of a apbara ia tba rarafiad flow xagina. Tha a^aatioaal baaia 
for tbia calcalatioa ia givaa ia Tabla 5,3. Tbia aabroatiaa calla DWNSTM to gat 
tba poat aoraal abock taaiparatare. TD, aad coapreaalbil ity factor. 2. 
Sabroatiaa HANSEN ia eallad to obtaia tba wiaooaity, XMU, at poat-abock praasara 
•ad tettp^rature lavala« 


Sabroatiaa STHEAT ia eallad by FAIRID. Tba argaaaata of tba CALL atataaaat 


ara: 


UI “ (iapat) Fraa atraaa waloeity 

RI " (iapat) Fraa atraaa daaaity 

T *■ (iapat) Fraa atraaa taaparatara 

P “ (iapat) Fraa atraaa praaaara 

AM “ (iapat) Fraa atraaa Maeb aaabar 

GANNA “ (iapat) Fraa atraaa apacifie baat ratio 

HI - (iapat) Fraa atraaa atatie aatbalpy 

TW ~ (iapat) Wall taaparatara 

Hf - (iapat) fall aatbalpy 

R - (iapat) Body radiaa 

QD - (oatpat) Coavactiwa baatiag rata 

B ~ (oatpat) Baat traaafar flla eoaffioiaat baaad oa taaparatara 

HTIL ~ (oatpat) Heat traaafar flla eoaffioiaat baaad oa aatbalpy 

HAW - (iapat) Laaiaar adiabatic wall aatbalpy 

TAW ~ (oatpat) Laaiaar adiabatic wall taaparatara 








AK2 - (ootpvt) Eftr«ftetioii paruaetts 

CBT “ (output) Stuntoa auabuu 

Lubulud eoMoau /FSSTN/ uad /DNSISN/ uro uaud to truaufot laforautloa foe 
tbo sboek ouleulutloa* 

5.4 SDBtOOTlNBS PMBBAT AND EBP 

Subroutiao PMHEAT la uaad to eoaputa fzoa Boloeular flow boatiag to aap 
vladward faeiag auzfaeo. Tbo fuaetioa program ERF ia tba orror fuaetloa call ia 
aubroutlao PMHEAT* Tha aquatioaa prograamad ia aubroutiaa PMHEAT aad fuaotloa 
EBP ara givaa la Tabla 5.4* 


Sobroutlaa PMHEAT la aallad by FAIBID* SPCHI aad SfCXL2. Tba argumaata of 

tba oall atataaaat ara: 

UI - (laput) Praa atraaa valoelty 

BI “ (laput) Praa atraaa daaalty 

T *- (laput) Praa atraaa taaparatura 

P - (laput) Praa atraaa praaaura 

AM - (laput) Praa atraaa Maeb auabar 

GAMMA ~ (laput) Praa atraaa apaolflo baat ratio 

El - (laput) Praa atraaa atatle aatbalpy 

TW - (laput) fall taaparatura 

BW - (laput) fall aatbalpy 

QC (output) Coaraotlva baatlag rata 

H - Output) Baat traaafar film ooafflolaat baaad oa taaparatura 

HTIL ~ (output) Boat traaafar film ooafflelaat baaad oa aatbalpy 

HAf ” (laput) Laalaar adlabatle wall aatbalpy 

TAf ” (output) Laalaar adlabatla wall taaparatura 

TH “ (laput) Tbata»6* aurfaoa aagla 

No ooaaoa block atataaaata ara uaad to traaafar laforaatloa* 

S.S 8DBBOOTINB SfClL 

Subroutlaa SfCIL la aallad by MAIN to ooaputa tba baat traaafar eoafflelaat 
aloag tba atagaatloa llaa of a awapt aylladar for laalaar or turbulaat flow* 
Tbla la doaa ualag aaplrleal ralatloaa baaad oa tba work of Cato for laalaar 
flow aad Jobaaoa for turbulaat flow aa glwaa la Tabla 5*5* Botb taebalquaa ara 
baaad oa aaplrlaal adjuataaata to tba atagaatloa poiat baat traaafar ooafflelaat 


U. ..«.. .» «. .,H.,„. „ ^ 

«....., ,„4i... U ^ ^ 

'»*« “•• p,„„. 

Mt« 1^ ».y„J4. .„,.^ g . ,,!/» J, 

...4 ».m „.„„ .. ^ 

. . .,U.4„ „4.„. «, .,. „.., ..,„, 

«.«. .M HHCX 44. 4„4.,.„ „.„ g... 

»«»..... w . 1 ^ n, „. ,^„. ,„g 
5.« SOBRODTINE SWCIL2 

1T4T'“ ““ " “ "•'• “‘* - - -. 

..4ifi.4 .0 .„.,..g g.,, ...ffl.,„„. 


An initial cheek ie aade on NONCON » miHA<\ ^ ^ 

on noncon - W(645) to dotermlao if r.rofi.d flow 

*• to be oonsldered. Subsequently, if noncom s a *v 

NONCON > 0, the subroutine stEO^xs Is 

..n.4 4P „c«.. „,.. gg ggg^ ^ 

UP... gg,. 

»..«...„ 1. .». ,g ,g. ^ 

^O. r |*®P»t) blinder rediut (ft) 

2 . T -««u*ir- “« - 

(...,.., f;:j|;:ran4„^.u.j...,. .„4.,.« 4... 


ii. m ■ I »o..4... i^ggg 

A 4.U4..P 1.,., M« ..iMUtio, „ „4, „ g„ 
4.44 4.44 f« „4 gj^^ ...gg^gg, 

Ur S.4P.4 pup...., 8.4roo4U. IIOLm i. t.jj,, „ 


otlbed in Table 5.6a and 
Per is set equal to 0,0 
obtain wall propertiea. 
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HANSEN !• to obtain woll Tlaooalty. FAYRID it collod to obtain tha 

apboc# ataination point boat tsanafar ooaffioiant. 


If IM - a Raaafiid FLa 

A raxafiad flow calanlation ia nada by ealling CYT. PCT ia aat aqnal to 

3,0 for ontpnt pnwpoaaa. 

If IRB » 3 Moleenlax Plow 

A ftaa aolaenlar flow oalonlation ia nada for a plata at tba local awaap 
angla. PCT ia aat aqnal to 10.0 for ontpnt pnrpoaaa. 


Tba labelad COMMCW baa baan addad to SfCfL2 and it ia eallad /FREMO/ • Tba 
wariablaa fron thia anbrontina tranaferrad tbrongb tba COMMCm/FREliO/ arc: 

AI2CIL •> Rarafaetion paraaiatar, 

XI ~ Stanton Nnmbar paranatar,^ 

plna tba wariablaa dafinad in aaetion 5.1. 

3.7 8DBR0DTINE CYT 

Snbrontina CYT ia naad to oonpnta rarafiad flow atagnation lino boating 
rata a for a right oironlar cylindar at arbitrary yaw anglaa. Tba aqnationa pro- 
graauDod in tbia anbrontina ara giwan in Tabla 5.7. 


Snbrontina CYT ia oallad by SWCYL2. Tba argnaanta of tba CALL 


ara: 


01 - (inpnt) 

RX - (inpnt) 

T - (inpnt) 

P - (inpnt) 

AM - (inpnt) 

GAMMA - (inpnt) 

HI - (inpnt) 

TW - (inpnt) 

HI - (inpnt) 

I - (inpnt) 

QC - (ontpnt) 

H Ontpnt) 

BTIL - (ontpnt) 


Fraa atraaa walocity 
Praa atraan danaity 

atraaa tanparatnra 
atraaa praaanra 

Fraa atraaa Macb nnnbar 

Fraa atraaa apaeifie boat ratio 

Fraa atraan atatie antbalpy 

Wall tanparatnra 

Wall antbalpy 

Diatanoa along tba tnrfaea 

Convactiwa boating rata 

Boat tranafar filn ooafficiant baaad on tanparatnra 
Boat tranafar filn eoaffieiant baaad on antbalpy 


Fraa 

Fraa 


\ 
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SAW - (input) LMiintt adiabatic wail entbalpy 

TAW - (output) Laainar adiabatic wall taaperatuca 

DLANB *> (input) Swoop anglo 

AK2CIL- (output) Rarof action Paraaetor* 

XI - (output) Stanton Nnabor Paraaetor, ^ 

Subroutine CYT oollo DWliSTM to oaleulato post noraal shook toaperaturos and 
prossuroa, Subsoquontly oubroutino HANSEN is oil led to ooaputo viscosity* 

Labeled eoisaons /FRSIli/ and /DNS1RM/ are used to transfer inforaation for 
the shock calculations* 

3*8 SDBROUTINB S9CYLi 

Subroutine SWC1L3 is also called by MAIN and eaploys the technique for 
coaputing the heat transfer coefficient along the stagnation line of an infinite 
swept oylinder as described in Table 5*8* This routine calls subroutines HOLIER 
and BINIRP* Hie arguaent of the CALL ststeaent for SWCYL3 contains the three 
following paraaeters: 

RN - (input) Cylinder radius 

ENGL - (output) Laainar heat transfer coefficient along the stagnation 
line 

ENCT — (output) Turbulent heat transfer coef'^icient along the 
stagnation line of an infinite swept cylindei* 

All other input inforaation necessary for coaputation is btou^ht into SWCTL3 

through ooaaon blooks* 

5*9 SDBROUTINB DEIRAL 

Subroutine DEIRAL is called by MAIN and is used in ealeulating the laainar 
and turbulent heat transfer coefficients about heaispherieal nose shapes* This 
routine enables the evaluation of circuaferential as well as stagnation point 
heat transfer coefficients* The laainar equations of Less* given in Table S*9a* 
are used to ooapute the laainar heat transfer coefficient* The turbulent equa- 
tions of Detra and Hidalgo, given in Table S.9b, are used to compute the turbu- 
lent heat transfer coefficient. 
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The arttMAntt of tlia CALL at 


for DE1RAL oostalo tba following |tar 


atara: 

IN - (input;) Sphara radina 

EL - (input) Running langth fron atagnation point 

FBI - (inpnt) Body angla (M dagraaa at tha atagnation point) 

ENGL - (ontpnt) Laainar Boat tranafar eeaffieiant 
BNCT - (ontpnt) Torbnlant boat tranafar ooaffioiant 

DETKAL aalla PAIBID to obtain tba laainar atagnation point boat tranafar ooaffi- 

aiant and aalla BINIltP to obtain tba diaaoeiation antbalpy, HD. 

5.10 SDBRODTINB BGOJtT 

Bnbrontina BOBBT ia aallad by MAIN and ia naad to calanlata tba laainar 
and tnrbnlant baat tranafar ooefficienta for flow owar a flat plata with nagli** 
gibla praaanra gradiant. Iba CALL atataaant for anbrontina BCKERT oontaina 6 


paraaatara: 


BLL 
BLT - 

ENL - 

ENT - 

ENCL - 
ENCT - 


(inpnt) Laainar running length (alraady eorraetad for 
croaaf low if eroaaflow option waa oboaan) 

(input) Turbulant running langtb (alraady oorractad for 
eroaaflow and virtual origin adjuatnant if tbaaa 
optiona wara oboaan) 

(input) Laainar Hanglar tranaforaatlon faotor to aeoount 
for axiayaaatrio flow eorraetion to flat**plata 
boating 

(input) Turbulant Hanglar tranaforaatlon factor to aoeount 
for aziayaaMtrie flow eorraetion to flat'^plata 
beating 

(output) Laainar baat tranafar eoaffieiant baaad on antbalpy 
(output) Turbulant baat tranafar eoaffieiant baaad on antbalpy 


Subroutine BXNIRP ia eallad to dateraina tba Prandtl nuabar PTj^ aa a funo- 
tion of tba Bekart rafaranea taaparatnra for laainar flow (raeovary factor la 
0.85) and tba ratio of tba boundary-layer edge praaanra to atandard ataoapbarie 
praaaura. Tba Bekart baat tranafar eoaffieiant for laainar flow ia ealenlatad 
uaing tba ralationa ia Table 5.10a. Naat anbrontina BXNTOP ia eallad to datar- 
aina tba Prandtl nuabar Pr^ aa a function of tba Bekart rafaranea taaparatnra 
for turbulant flow (raeovary faotor ia 0.88) and the ratio of tba boundary-layer 


pr««sor« to ttandard •taosphoric prtasnr*. Tha torbnUnt flow baat 
tranafor coafflolant ia oalcalatad aalng tba ralatioaa glvaa in Tabla 5.10b. 


If t:ao paraaatar NERROR ^ 0 la raturnad Iron BINIRP, than ECKERT rat *rna to 
MAIN and prooaada to raad ia tha aaxt data caaa. Ilia Ecfcart rafaranca propar- 
tlaa Uaainar aad tarbalant) axa oalcnlatad la aabroutina FLOW baaad oa a con- 
ataat Pcaadtl anabar (laalaar raaovary factor la 0.85 and tba tarbalaat lacovary 
factor la 0.88). Oalnf tbla conataat Fraadtl naabar tba Eckart rafaraace taa- 
paratuza la dataralnad. Tbla rafaranca taaparatnra la than aaad la conjaactloa 
with BINIRP to dataralna tba Prandtl naabar (now coaaldarad aa a function of 
taaparatnra and praaauro) wblcb waa originally taken to ba a conatant. Strictly 
apaaklag, tba Prandtl nnabar abould ba dataralnad by an Itaratlwa procadora. 
Rowawar. tbla probably does not cauaa any appraciabla error. 

S.U SDBROOTINES SPCSI AND FSUBC 


Subrontlna SPCBI la called by MAIN aad la aaad to calcalata tba tarbalant 
heat tranafar coefficient for turbulent flow, and alao tha Eckart rafaranca 
taebni^aa la need to coapata tba laminar beat tranafar coefficient. Bot!, cor- 
raapond to flow over a flat plate with a negligible praaaara gradient. Hie CALL 


for aabroutina SPCBI contains aavan paraaatara. 


ELL 

ELT 


ENL 


ENT 


ENa - 
ENCT - 
RANPLO - 


(input) Laminar running length (already corrected for 
crossflow If crossflow option was cbosan) 

(input) Tarbalant rannlng length (already corrected for 
crossflow and wlrtaal origin adjastmaat If tbaaa 
options ware chosen) 

(Inpat) Laminar Mangier transformation factor to aoooant 
for aslsynmatrio flow corractlon to flat-plata 
beating 

(Inpat) Tarbalant Mangier transformation factor to account 
for anisymmatrle flow correction to flat-plata 
beating 

(output) Laminar beat transfer coefficient baaad on enthalpy 
(ontpat) Turbulent beat transfer coafflclaat baaad on enthalpy 
(Inpat) If a walaa greater than aero la Input into location 
919, than tba Von Karmsn fora of tba Reynolds 
analogy factor la used. 
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As isltial ohaek is sad* os NONOON W(64<S) to datanaisa if bosodarj lajar 
flow ia to ba aooaidarad (jaa If NON CON > O), Sobaa^saotljf it NONOON > 0« tha 
esbaontisa BBQUIB 1* oaltad bj SPCHl, Baaad os tba walsa of ItB* tba ostpst 
from RBOIKB, a boosdarp layar« rarafiad oz fra* solaoulai )ala«latios It sada* 

II MB ■ 1 BOMtfUT Ufft rigf 

A boQodazj lapar flow calcalatioo it sada a* daacribad la Tabla S.ll for 
tsrbolast flow asd Tabla 5.10a for laslaar flow. PCT it tat aqsal to 0.0 for 
ontpnt parpotta. 

SPC8I eallt asbrostisa FSUBC to eospnta F for a raal gaa saint BOLIEB. 

Tba arfsaast of FSUBC oontaina flva parasetart: 

HE - (ispst) Batbalpp at tba tdga of tba boundary layar 
BW “ (input) Wall antbalpy 

HKECT - (input) Turbulant Eoktrt raftranot antbalpy 
PE - (input) Boundary**layar adga praatnra 
FCINV - (output) 1/F 

0 

If no problas ia anoonntarad wltb tba Nollier diagras. than F it rtturnad to 

c 

8PCBI. If tba tbarBOdynasle pzopoztiaa oannct ba dottrninad by MOLIEB* than 

IDEAL • 1 ia raturnad to SPCSI and an idaal gat ealenlatlon of F it undartakan. 

0 

Having obtainad F^» tba turbulant baat tranafaa co«ffioiant la than datarsinad 

froB tba Spalding-’Cbi equations wbieb bava baan ebeekad and found to ba corraat. 

If a valua graatar tban aaro ia atorad in location 319# tban tba Von Eaman Bay- 

noldt analogy factor ia uaad in datarmining tba turbulant baat trantfar cocffi** 

ciant. Bowavar* if taro ia atorad in location 919. tban tba Eaynolda analogy 

2/3 

factor ia uaad at Pr • 

Aftar tba turbulant baat trantfar coafficiant bat baan datarsinad* 8PCHI 
naat procaada to coaputa tha laslaar baat trantfar coafficiant by tba Eckart 
rafaranca tacbniqua. Both tba iasinar and turbulent baat trantfar coafficiantt 
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•r« eorr«ot«d for cxottflow and for •sltynaotrio flow offoeto (via tha laainar 
aod torbalaat Maaglar traaaf oraatloa factors) bafora balaf aaad la tba boat 
traaafor ooafflelaat aqaatloaa. 

If lEE ■ 2 laraflad glow 

NFCS > 0.0, NPC8 - W(314) 

If NFCS > 0.0, thoa SPCBI oalla SPFP. 

Tba raraflad flow boat traaafor eooffloloat for a abarp flat plata at 
raro aaglo of attack la ooapated, PCT la act eqaal to 4.0 for oatpat 
parpoaea. 

HPCS 1 0.0 

If NFCS i 0.0 tbea SPCDI calls CONE. Tba rarefied flow beat traaaf er 
eoefflcleat for a abarp aose eone la ooapated. IfCt la set eqoal to 5.0 
for ootpnt parpoaea. 

IIIEE « 3 Free Molecalar Plow 

A free aoleealar flow beating oaloalatloa la aade for a plate at tbe loeal 
bodjr aagle by oalllng FMHEAT. PCT la aet eqaal to 10.0 for oatpat parpoaea. 

5.12 SUBROUTINE SPFP 


Sabroatlae SPFP la aaed to ooapate flat plate beating rates In tbe rarefied 
flow reglae. Tbe eqaatlons prograaaed In tbla aabroatlne are given In Table 
5.12. 


Sabroatlae SPFP la called by SPCBI. Tbe argaaeats of tbe CALL atateaaat 


are: 


UI 

- 

(lapat) 

Free 

R! 

- 

( lapat ) 

Free 

T 


(Input) 

Free 

P 

- 

(lapat) 

Free 

AN 

- 

(Input) 

Free 

GAMMA 

- 

( Input ) 

Free 

BI 

-- 

(Input) 

Free 

TW 

- 

( lapat ) 

fall 

BW 

- 

(Input) 

Wall 


atreaa veloelty 

atreaa denalty 

atreaa teaperatare 

atreaa preaaare 

atreaa Maeb aaaber 

atreaa apeelfle beat ratio 

atraaa atatle eatbalpy 

teaperatare 

entbalpy 


X “ (input) Pittance along the tnrfaee 

QC ” (outpnt) ConvcetiTC beating rata 

H - (ontpnt) Heat tranafer fila coaflieiant baaed on taaparatnra 

HTIL - (ontpnt) Beat tranafer fila eoaffieiant baaed on enthalpp 

HAH *■ (input) Laainar adiabatic vail entbalpp 

TAW “ (ontpnt) Laainar adiabatic vail teap eratnre 

<SMI - (ontp.nt.).Aaxnfantion.paraaeter, H /C^/Re 


Subroutine 8PFP calla HANSEN to coapnte viacoaity for nae in calculating 
tbe CbapBan''Bubeain conatant. 

No coaaon atateaenta are uaed« 

S.13 SDBBODTINE (XPIE 


Subroutine CONE ia need to coapnte rarefied flov beating rate a for a abarp 
cone. The eqnationa prograaaed in tbia anbroutine are given In Table 5,13, 

Subroutine CONE ia called by SPCHI, The argument a of tbe CALL atateaent 


are: 

UI - (input) 
HI - ( input ) 
T - (input) 
P ~ (input) 
AM - (input) 
6AMIA - (input) 
HI - ( input ) 
TW - Unput) 
HW - (input) 
X ~ (input) 
QC - (output) 
H ~ (output) 
HTIL (output) 

HAW “ (input) 
TAW - (output) 
TCP (input) 

CHI (output) 

ZBTA “ (output ) 
N8B “ ( input ) 


Free atreaa velocity 

Free atreaa denaity 

Free atreaa teaperatnre 

Free atreaa preaaure 

Free atreaa Na oh nuaber 

Free atreaa apeoific heat ratio 

Free atreaa atatic enthalpy 

Wall teaperatnre 

Wall enthalpy 

Piatance along tbe aurface 

(k>nveotive beating rate 

Heat tranafer fila coefficient baaed on teaperature 
Heat tranafer fila coefficient baaed on enthalpy 
Laainar adiabatic vail enthalpy 
Laainar adiabatic vail teaperature 
Seaivertex angle of tbe cone 
tarefaction paraaeter, ^ 

Stanton Nuaber paraaeter, C 
NSB - W(6S0) 

NSB “ 0 Sharp cone curve fit 
NSB M 1 Blunt cone curve fit 
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SvbroBtla* CONE eallt OWNSTM to etlcolato post sormal shock tcaperstnrc and 
prasaora. SnhaaqBa&tlj aobroBtiaa HANSEN is eallad to eoapata Tlacoaity basad 
OB a rafaranca taaparatnra. 

Labalad coanona /FRSTM/ aad /DNSIBM/ ara naad to tranafax infoxnatioa for 
tha shock calcalatloa, 

5.14 SOBRODTINB SBOWJR 

Sabroatiaa RHOMUR is called by MAIN aad is ased to eoapata laaiaar and tar“ 
balant hast tranafar coaffieiaats ovar a flat plate aaiag tha ralatioaa in Table 
5.14. Sabroatiaa BINIRP is called by RHOMUR to ooapata tha Pzandtl naabar as a 
faactioa of rafaranca taaparatara and as a fanction of tha ratio of tha 
boandazy-layaz edge pzassaze to standard atmosphazic prassara. Also BINTRP is 
eallad to ooapata tha specific dissociation enthalpy as a fanction of zafazaaca 
taaparatara aad as a fanction of tha ratio of tha zafazaaca density to standard 
ataosphazie density. Tha azpressioas for tha laainaz and tazbalant heat 
transfer eoafficiant haaa bean soaewhat siaplifiad fzoa tha original •*- 

pzassions. Tha aqaivalant zanning length is datarainad fzoa CRSFLW ar.a tha ozi- 
ginal crossflow azpzassions are not asad. Also there Is nc adjastaant to 
tha heat transfer eoafficiant to accoant for the inezaasad heating daa to a*- 
isynaatzie flow (Mangier tzansfozaation factor). Tha argar:^nt of tha CALL 
atataaaat for RHOMUR eontaias tha foar following paraaatara: 

~ (inpat) Laaiinar zanning length which nay have bean eorractad 
for crossflow affects 

~ (iopnt) Tarbalaat zanning length which nay have bean 

eorractad for crossflow affects and wirtaal origin 
adjastaant 

ENO< ~ (oatpat) Laaiinar heat transfer eoafficiant 

ENCT — (oatpat) Tarbalaat hast transfer eoafficiant 

All other inforaation necessary for eoapatatioa is broaght into RHOMUR throagh 


coaaon blocks 


I 

I 5. IS 8UBR0DT1NB LE81A8 

j SttbroQtiM LB81DE It called froa MAIN after a call to FA^ID* to ealeolate 

the averaae heat traaefer eoefflelent to the leeeide of an arbiter. The arerafe 

heat transfer eoeffleieat is ealealated aeiat the relatione in Table 5.15. This 

sabrontiae ie called If NHFLAG ■ 9, The arfOBenta of the call are: 

IN - (iapat) Radiae of scale, fall scale RN • 1.0 ft. 

BVD - (lapat) W(21). Windward fall enthalpy (Bta/lbn) 

BN(X “ (oatpat) Averape tarbaleat or Isaiaar leeward side heat 
transfer coefficient 
ENCT - (oatpat) ■ ENCL 

The cosmon FLWFLD ie ased to transfer post aoraal shock aalaes obtained by the 
call to FAIRID. 

5.16 SUBROUTINES FLAFH, DEFL, DSTNL. CPIAF AND BTHICX 

Sabroatine FLAPH is oalle fr<« MAIN when NHFLAG - 10 to calcalate the peak 

reattacbaent heating to a flap. The eqaations ased for this calcalatioa are 

glaen in Table 5.16. The argaments of the call are as follows: 

EL - (inpat) Rnnaing length 

FLENO - (iapat) Flap length 

FANG > (iapat) Flap angle 

TV ~ (inpat) Wall teaperatara 

BT > (iapat) Total enthalpy 

ENCLO - (oatpat) Laminar heat transfer coefficient 
ENCID “ (oatpat) Tarbalent heat transfer coefficient I 

IQUIT - (oatpat) Flag IQUIT ■ 0 Separation occars 

IQUIT • 1 No aeparation occars 
PCI - (oatpat) Transition percentage 

A cMBon etateaent. FLPETB. ia ased to transfer the following edge properties 

froa NAIN into the sabroatine: 

EDOEP(l.I) > Preseare 
E0G£P(2*X) - Teaperatare 
EDOEP(d.I) - Density 
ED0EP(4»I) “ Speed of eoand 
ED0EP(5.X) - Velocity 
ED6EP(6«I) “ Mach naaber 
ED0EP(7.I) - Viscosity 
EDOEP(S.l) - Specific heat ratio 

where Z • 1 - Before the flap shock edge conditions 
I 2 * After the flap shock edge eonditiona 


i 
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Sobrootlfio FLAPH e«ll« BIBICK to obtiln a laalaar and torbalent boundary 
layor tbieknoaa* This oaienlatioa it aadt ntlng Ecktrt rtfortnet oondltlont in 
tbo followiny oqnationa: 

6 ■ 5.2X/»4#^* t Laainar boundary layer thieknett 

® “ 0,1541/ (Ro^ ) ' t Turbulent boundary layer tbiekneaa 
Next aubroutioe CPIAF ia oalled to oaleulate tbe incipient aeparation preaaure 
ooeffieient, lltia ia followed by a call to OEFL to determine if tbe input flap 
angle ia aufficiently large enough to produce aeparation. If the input angle la 
inaufficicnt to cauae aeparation, IQUIT ia aet equal to 1 and the aubroutine re- 
tnrna to the MAIN, If the aubroutine returna to the MAIN with IQUIT • 1, then 
the heating ia calculated by a call to SPCBI, 

If aeparation ia found to occur, the plateau preaaure ia calculated and 
OEFL ia called to calculate the dividing atreamline angle. Next the aeparation 
length ia calculated by a call to OSTHL, If the aeparation geometry indieatea 
that impingement will not occur for the input flap length, IQUIT • 1 and the av 
broutine returna to the MAIN. 

If impingement occura, the ehear layer thickneaa ia calculated. Subroutine 
BANSEN ia called to calculate the wall viaeoaity. Finally, the heat tranafor 
coefficient ia calculated for eithor laminar or turbulent flow with PCT aet 
equal to 0,0 or 1,0, The laminar and turbulent heat tranafer coefficicnta are 
equated ainee thia routine can calculate only one type for a given flowfield. 
Note that in MAIN if NBFLAO • 10 and IQUIT • 0 then the call to TRANS ia 
akipped, 

5,17 SUBROUTINES FINH, FINPKH AND FINPKP 

Subroutine FINB ie called by MAIN when NHFLAG “ 11 to calculate tbe peak 
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h«ttiat prodaetd by a fin oa tba adjaeaat tarfaea* Tba ataatioaa 
aaad for tbia ealealatiea aro giroa ia Tabla 5.17 Tba argaaaata of tba oall 
ara aa followat 


ALP - 
XP 
IP 
BL 
IW 
10 

AMPPT - 
AMPPL - 
AMniL - 
ANFIT - 
BNCLO - 
ENCID - 


(iapat) 

(lapat) 

(oatpat) 

(lapat) 

(iapat) 

Uapat) 

(oatpat) 

(oatpat) 

(oatpat) 

(oatpat) 

(oatpat /iapat) 

(oatpat/iapat) 


Bffaatiao fia aagla of attaak 

Laagtb aloag fia sarfaoa 

Noraal distaaoa item fia to peak baatiag 

Baaaiag laagtb to fia laadiag odga 

Wall taaparatara 

Total aatbalpy 

Tarbaleat praatara aaolifioatioa 
Laaiaar praaaara aa^ fiaatioa 
Laaiaar baatiag aaplafiaatioa 
Tarbalaat baatiag aaplifiaatioa 
Laaiaar beat traaafar aoaffieiaat 
Tarbalaat boat traaafar aaaffiaiaat 


Tba aoaaoa atataaaat. FLPHIG. ia aaad to traaafar adga proportiaa apatraaa of 
tba fia aboek froa MAIN iato tba aabroatiaa. Iba aariablaa trauafarrad a.ra da~^ 


fiaad ia aabaaatioa 5.16. 


Sabroatiaa FINE ealla BTHICK to obtaia tba laaiaar aad tarbalaat boaadary 
layar tbiekaaaa aaiag tba ralatioaa givaa ia aabaaatioa 5.16. Mart tba aboak 
aagla prodaead by tba fia ia eonpatad aaiag a aall to PCSW. If tba aboek ia at- 
taabad. tba paak praaaara aaiplif ieatioa ia aai6ul*>tad for laaiiaar aad tarbalaat 
flow aaiag two ealla to FINPKP. Tba paak baatiag aaplifiaatioa for laaiaar aad 
tarbalaat ia tbaa aoapatad aaiag two ealla to FINPKH. Tba iapat baat traaafar 
aoaffiaiaata ara aodlfiad aaiag tba baatiag aaplifiaatioa faatora. 


Sabroatiaa 8PG1I ia alwaya eallad bafora tba aall to FINH to obtaia tba aa** 
diatarbad baat traaafar coaffloiaat. Nota that ia MAIN, if NHFLAO •* 11. tba 
baatiag aaplifiaatioa faetor aad adga praaaara ara oalealatad baaad oa PCX after 
tba aall to TBANS. 

5.18 8WBR0DT1NB BADEQT 

Sabroatiaa BAlffiQT ia eallad by MAIN aad ia aaad to aoapata tba radiatioa 
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noillbrin teBp«r«tura. Hie radiation aqoilibriun tanparatnra ia datarninad 
via tha Nawton-Raphaon itaration taabaiqoa. Tba wall tanparatnra ia naad a a an 
initial maaa and eonyarianea to tba radiation aqnilibrinn tanparatnra la conai- 
darad to ba aabiavad aftar 50 itarationa. or aftar two anccaa.iwa itaration. ara 

within 0.5 parent of ona anothar. Tha first of that, two condition, to oeenr 
dafinaa eonvarganea. 


Tha argonant of tha CALL atatanant for anbroutina RADEOT contain, tha fol- 
lowing fiva paranatara: 


ENC 

HR 

EHIS 


(input) 

(input) 

(input) 


W - (input) 


TRE - (output) 


tranafar coaffieiant 
Raeovary anthalpy 

Effactiva anisaivity - product of ahapa factor and 
aniaaiwity 

Conduction or thin-ahin wall tanparatura which is 
uaad as an initial guess to tha radiation 
aquillbriun tanparatura in tha iteration achana 
Radiation aquilibriun tanparatnra 
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Table 5.1 

STAGNATION POINT HEATING 


OhSuiiViAL 

OF POOR QUALITY 


Stagnation point heating Is based on the theory of Fay and Riddell. 




where 


!!u 


0.76 Pr^j /ptnA I' l*-e ■ ') Si 


/^\ • 1 JOj (Pj- O ^0'^ • 

^arj ot 

■ 1 12a ( P • P ) for a cylinder and M^>1.12 

TrV C.1^ 


Pt 

U (1.5 - 0.378M* - 0.02625 ) for a sphere and M^< 1.22 

* 00 00 

R 

(2.0 - 0.872M* - 0.328MJ^ ) for a cylinder and »1.< 1*12 


and the dissociation enthalpy 

«d • ^o^’o^W 

hj ■ 6636.26 Btu/Lbm 

h* ■ 14456.53 !'tu/Lbm 
N 


Nomenclature 

q • 

0 * 

^c 

9 ■ 

P • 

V ■ 

u • 

p • 

H ■ 


Heating Rate (Btu/ft*sec.) 
1 (slug ft./lb^sec.*) 

32.174 (Lbm/slug) 

Density (s1u9/ft,*) 
Viscosity (Lbxsec/ft*) 
Velocity (ft.;sec.) 
Pressure (Lb^/ft.*) 
Enthalpy (Btu/lbm) 
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Velocity Gradient 


Table 5.1 (Cont. 1) 

ORIGINAL RAGS Pi 
OF POOR QUALITY 

Lg ■ Lewis Number 

• Dissociation enthalpy of air (Btu/Lbm) 

• Mass fraction of species i (i • 0, oxygen and i » N, nitrogen) 

R ■ Nose radius (ft.) 

Pj, ■ Prandtl Number 

0 ■ Sphere Body Angle 

Subscri cs 

e » Edge 

t ■ Stagnation (post shock) 

w « Wall 

• » Free stream 




O 

-D 






References 


Fay, J. A. and Riddell. F. R., "Theory of Stagnation Point Heat Transfer 
lo Sci. . Vol. 25, No. 2, Feb. 1958, pp. 

7o*o5 f 1 21 » 



Table 5.2 

BOUNDARY LAYER, RAREFIED AND 
FREE MOLECULAR FLIGHT REGIME CRITERION 


A - — = 


Flight Regime Selection Parameter 


If 

A < 

0.05 

Boundary Layer 

If 

0.05 

< A < 3.0 

Rarefied 

If 

A > 

3.0 

Free Molecular 


where 

■ Free Stream Mach Number 

Re^ • Free Stream Reynolds Number Based on Running 
Length or Radius 

2e ■ Post Normal Shock Compressibility 


Table 5.3 


ORIQIWAL PAGS 
OF POOR QUALITY 


RAREFIED FLOW STAGNATION POINT HEAT TRANSFER EQUATIONS 


The rarefied flow heating to the sphere stagnation point based on the work 
of Engel and PraharaJ Is as follows. 


(1) % • ■** (Reference temperature) 


U) T,.T 




(Free stream stagnation temperature) 


(3) K* ■ e (Rarefaction parameter) 



where e ■ 


Id 

2Y 


(4) Heat transfer coefficient 

2 


1-0 

a > -0.235256 

0 

a^ • -0.303095 

a ■ -0.0779538 
a 

(5) Heat transfer 
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TabU 5.3 (Cent. 1) 


ORfQiMAL PAGE j;a 
OF POOR quality 



Nomenclature 




Subscripts 


I«K» - feP ^Tr 
Ur«»T« 

SUgiMtIon Point Hoot Tronsftr In Low Riynolds Nunbor Flow 

Stanton Number 
H Enthalpy 

Free Stream Mach Number 
T Temperature 

U Velocity 

2 Post Normal Shock C(»npress1b111ty 

P Density 

y Specific Heat Ratio 

V Viscosity 

• Free Stream 

W Wall 

6 • e Post Normal Shock 

0 Total 


Reference 


* S.C., "External Tank Rarefied Aerothermodynamics," 
REMTECH Inc. RTR 022-1, January 1978. 
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Table 6,4 (Cont. 1) 



ORfQ;^!'V, , 

C»F foou QUAL-ri 


( 6 ) ^ • j ^ 


I 



These equations are based on the work of Oppenhelm. 


Nomenclature 

k - Boltzman's Constant 
M« - Freestream Mach Number 
P - Freestream Pressure 
q • Heating Rate 
- Wall Temperature 
T» - Free Stream Temperature 
U« - Free Stream Velocity Heat Ratio 
Y - Free Stream Specific Heat Ratio 
a - Accommodation Coefficient 


References 


Abramowitz, M., andStegun* !. A,, eds» Handbook of Mathematical Functions. 
Dover Publications, New York, 1965. 


Oppenhelm, A. K., “Generalized Theory of Convective Heat Transfer 
Free-Molecule Flow", J. Aeron. Scl.. Jan. 1953, p. 49. 


In A 




Table 5.5 

CATO/JOHNSON SWEPT CYLINDER HEATING 


PriC. .; jrj 

OF POOR QUALITY 




For laminar flow, 


hj^ ■ 0 




For turbulent flow. 


hj • 0.75-^ h [ -0 ^u*^n ]^*^ I 

0 I I |0.01714 ♦ 0.1 

The «1ttp»c.t.on fetor. K^. , eluded t„ 


01235 sin ( 3.53 (A-10*: 


the definition of h . 

0 



A, Sweep Angle (Deg.) 



Table 5.5 (Cont. 1) 


Nomenclature 


Stagnation Point Heat Transfer Coefficient (Sphere) 
h|^ Laminar Heat Transfer Coefficient 
Turbulent Heat Transfer Coefficient 
K Multiplier Factor 

Cylinder Radius 
Vy Velocity (Upstream) 

Py Density (Upstream) 

A Sweep Angle 

Py Viscosity (Upstream) 


References 

Cato, 6.C., "Heat Transfer to the Leading Edge of a Yuwed Wing", 
Memorandum A*260>TH>57*115. 

Johnson, W.A., "Turbulent Heat Transfer to a Yawed Circular Cylinder", 
Manorandum A2-260-TH-59-218, September 1959. 


•0'C>XC:»'0T3 3-«0 


ORIGINAL PAGu ?'? 
OF POOR QUALITY 


Tftbis 5»6® 

BECKWITH AND GALLAGHER TURBULENT REAL GAS 
YAWED CYLINDER STAGNATION LINE HEATING 




Nomenclature 


32.174 Lbm/slug 

Heat Transfer Coefficient Lbm/frsec 

Pressure Lbf/ft 

Prandtl Number @ Edge Conditions Dimensionless 
Radius of Cylinder Feet 

Velocity Ft/sec 

Distance from Stagnation Line Feet 

Sweep Angle Degrees 

Viscosity Lbm/ft-sec 

Density Slugs/ft’ 


Subscripts; 

SL 


Stagnation Line 
Edge Conditions 



Evaluate at Eckert's 


Reference Enthalpy & Edge 
Pressure 


Total 



A, Snmp Angl* (Dtg) 


Reference; 

Beckwith, I. E., and Gallagher, J. J., "Local Heat Transfer 
Temperatures on a Yawed Cylinder at a Mach Number of 4.15 and High Rey- 
nolds Number, NASA TR R-104, 1961. 


-r. 



Table 5. 6b 

LAMINAR CYLINDER STAGNATION LINE HEATING 


h - 0.707 hj Cos' * A 


h, • Sphere Stegintlon Point Keat Trensfer Coefficient 

h • Cylinder Stegnetlon line Heat Transfer Coefficient 
A ■ Sweep Angle 


Laminar Flow 
Cylinder Data 
■3.7 From 

Brewer, MSFC 


A, Sweep Angle (Deg.) 



Table 5.7 


ORIQINAL P/'jC,'^ 

OF POOR QUmLI'iT 


RAREFIED FLOW YAWED CYLINDER 
STAGNATION LINE HEAT TRANSFER EQUATIONS 


The rarefied flow heating to a cylinder based on the work of Engel and 
Praharaj Is as follows: 


(1 ) ^ (^w * cos*A) 


(2) Re« ■ -BeHasL 

Vm 

(3) C*«iiLls 2 - 

VWTr 


Re. Ze 


(4) It* • 

* ♦ P.s1n*A/p u «) •• ; 


(5) C • C|^/(cos*A + P^s1n*A/p^U^*) •* 

(6) Iteat Transfer Correlation 


’•9.. ? ■ ,1, ai 009.^» )< 

■ - 0.377656 ■ 

• - 0.368580 
da * - 0.0461064 

(7) Heat Transfer 

<l • 0- Cu (Ho - H^) 
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Tab!® 5.7 (Cent. 1) 


vWQINAL P/teS 
OF POOR QUAUIY 





iiHu&afiuuiu 
^fft-aaBlkii - c^i -?r 

iKiriefbdh;^ 

sssssimii 

nmaiBi 


■•■maBMai! 
MffiSijEnpSiiinai 





, JBgiaSilgt!iBaBSil88gil 

iliiyiiiiiiiiiilliiiil^ 




A£0C DATA 

Kim Gtugt A Conf. » 


o 4635 
0 4636 
O 4637 
□ 4682 
Q 4863 


29 60* ET .177 

29 JO ET .177 

29 0 CT .177 

29 0 n .177 

S"C;yl .177 
. ? S"fyl .} 77 « 



Nomenclature 


Subscripts 


Reference 


Cylinder Stagnation Line Heat Transfer Data with Sweep Data Included 

I Stanton Number 

Enthalpy 

, Free Stream Mach Number 

Temperature 
Velocity 

Post Normal Shock Compressibility 
Density 

Specific Heat Ratio 
Viscosity 


Free Stream 
Wall 

Post Normal Shock 
Total 


e • edge 


wilTECT‘iM.*RTIl'’o2^lt^j.nM;; 1 I 78 ?™*' Aeroth.n.odyn«Mes 
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Table S.«} 

SWEPT CYLINDER HEATING 


OF POO« QljALi r/ 


The "p«vJ' fuethod for swept cylinders Is defined as follows: 
r r 


H^c.o ■ ♦ ««> 


^mc,o " ^^^mc,o* 


0 p 


mC(0 


r« * 0 
0 


.55 

•K^o) 


-0.5 


- 0.90 /'*St.9'‘st.9\ 

V »w'‘w / 


0.5 


f(P|^. T^) 


L,o 


[’ ♦ ■'o] ^ 


‘T.o " [ 


1 + 0.77 r, 


1 

>J r 


’ f 


2(Pj - P„) 
^ 


Xeq 


9 K T 


Xeq 


T 9O 51/ -r 

®'^r^T,( 


Xeq 


T.o 


0.333 


•"l.o 


'^stag^stag 


PcV 


ri 
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Table 5.8 (Cont. 1) 


r '■ 1 ^ 

DE POOr< QUtiLliY 


H»c ■ «»c.o * '>•'“‘<'’1 - 


Pmc ■ ^*"inc’ ’’i* 


•, . [o.96(e/* -O.s] 


\ • (1 ♦ r^)p/ 


(T, ♦ 2000) 

PrUfV^ML,^ 

?\ Wj« (Xeq^) 


0.370 


[logio(Re,.cyl * 


• “o,f • 


i . 1 ♦ (Le°*®^ - 1) 
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Table 5.8 (Cont. 2) 


ri;>SQiK*AL FA0.5! 

OF POOR QUALn> 


For laminar flow, 


h 


L 



Re 


r.cyl 


Cf 


r.L 


For turbulent flow. 


0. 5g 



Re 


r.cyl 


Cf 


r,t 


Nomenclature 

Cf 

h 

H 

Le 

P 

Pr 

Re 

? 

V 

P 

y 

Subscripts 

L 

r 

stag 

f 

t 


Skin Friction Coefficient 

Heat Transfer Coefficient 

Enthal py 

Lewis Number 

Pressure 

Prandtl Number 

Reynolds Number 

Cylinder Radius 

Temperature 

Velocity 

Density 

Viscosity 


Laminar 

Eckerts Reference 

Post Normal Shock Stagnation 

Stagnation Line 

Total 

Turbulent 


References 

Nagel et al., "Analysis of Hypersonic Pressure and Heat Transfer Tests 
on Delta Wings with Laminar and Turbulent Boundary Layers," NASA 
CR-535, August 1966. 


Thomas et al . , "Advanced Re-entry Systems Heat Transfer Manual for 
hypersonic Flight," Technical Report AFFDL-TR-65-195, October 1966. 
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Table 5.9a 


OF POOR QUAU ry 


LEES' HEMISPHERICAL DISTRIBUTION METHOD FOR LAMINAR FLOW 


Lees' heating distribution over a sphere for an Ideal gas Is expressed 


as: 


0(e) 


Nomenclature 




20s1ne 


— 7W 


co8*e 


* -m 


(’ ■ w) [ 

r 

rM«2 


0s1n40 + 1 - cos40 
2 §“ 


+ YM«2 I 0^ - 9s1n20 + 1 - cos20 

2 


0 


] 


h • Heat Transfer Coefficient 

hj 5 p*-« Stagnation Point Heat Transfer Coefficient 

M» ■ Free Stream Mach Number 
y » Free Stream Specific Heat Ratio 
0 ■ Local Body Angle 


Reference 

Lees, Lester, "Laminar Heat Transfer Over Blunt-Nose Bodies at 
Hypersonic Flight Speeds". Jet Propulsion. April 1956. 
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where 

* 1 *’ 

1 (6 < 2S”) A. (6 > 25*) 

1 1.03754 0,9645? 

2 0.0043776 0.01107 

3 -0.6187 X 10"^ 0.842558 x 10"^ 

Nomenclature 

Gravitational Constant 
Turbulent Heat Transfer Coefficient 
Enthal py 

Turbulent Multiplier Factor 
Running Length 
Prandtl Number 
Velocity 
Density 
Viscosity 

Sphere Tangency Angle (90 Degrees at the Stagnation Point) 

Dissociation 
Edge 
Total 


Detra, H.W. and Hidalgo, H., "Generalized Heat Transfer Formulae and 
Graphs", AVCO Research Report 72, March 1960. 
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0 


e 

T 

Reference 


Table S.lOa 


ECKERT'S LAMINAR FLAT PLATE HEATING METHOD 


Eckert's laminar flat plate heating relation can be expressed as: 

h • 0.332 9c D*u« ... 

Obm/ft*sec) 

where 

R* m Reynolds no. evaluated at 

e V* reference conditions 

The reference conditions are obtained from the reference enthalpy 
H* « He + 0.5 (litf - He) + 0.22(Ha„ - He) 
and edge pressure. Other reference properties are evaluated as: 
p* » f (H*. Pe) 
u* * f (H*. Pe) 


For geometries other than a flat plate, the Mangier transformation 
may be evaluated as: 

% • 1 Flat Plate 

■ 3 Cone 

Nomenclature 

g^ » 32.174 Ibm/slug 

h ■ Heat Transfer Coefficient (lbm/ft*sec) 

He » Edge Enthalpy (Btu/lbm) 

H-,. « Adiabatic Wall Enthalpy (Btu/lbm) 

• Wall Enthalpy (Btu/lbm) 

Pe = Edge Pressure (atm. ) 

Ue “ Edge Velocity (ft/sec) 

X • Surface Distance From Origin to Point of Interest 
p* ■ Reference Density 
M* • Reference Viscosity 
Tn, “ Mangier Transformation 



Nomenclature 

9c ■ 32.174 Ibm/slug 

h - Heat Transfer Coefficient (Ibm/ft^sec) 

Pr* ■ Prandtl Number at Reference Conditions 

Ug » Edge Velocity (ft/sec) 

X - Surface Distance From Origin to Point of Interest 
P* ■ Reference Density 

P* ■ Reference Viscosity 

t|n * Mangier Transformation 


Reference 






Table 5.11 


SPALDING-CHl MnHOD MODIFIED FOR REAL UAS HEAT TRANSFER 
FOR TURBULENT BOUNDARY LAVER FLOW 


(1) Spfild1ng-Ch1 define F^, Fj^, Fjj^ (Spalding and Chi) which are 
functions of Mach Number and Temperature clone such that 
1/2 C^Fg ■ '^e(’^RQReg) " V'sCFr^ROj) 


(2) F^ 







Where g(l) » 9.2808635 
g(2) • -4.7340248 
g(3) - 6.6858663 10"* 
g(4) ■ -4.1876614 10‘* 
g(5)'- -5.5054577 lO"* 
g(6) • 2.8367291 10“*» 
9(7) - -2.1249608 10“ • 
9(8) • 8.0162000 10“’ 
g(9) • -1.5900985 j0“* 
9(10) - 1.3236350 *10“ “ 



■ Mangier Transformation (Komar) 
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Nomenclature 


Table 5.11 (Cent. 1) 


^TURB 

H 

P 

Pr 

Reg, Re 


s 


S 

St 

T 

u 

Z 

P 


Skin Friction Coefficient 

Heat Transfer Coefficient For Turbulent Flow 

Enthalpy 

Pressure 

Prandtl Number 

^ocal Reynolds Number Based on Momentum 
Thickness 0 and Characteristic Length s 
Respectively 

Inverse of Von Karman Reynolds Analogy Factor 

Stanton Number 

Temperature 

Velocity 

Compressibility 

Density 


Subscripts 

aw 

e 

1 

w 


Recovery 

Local 

Incompressible 

Wall 


References 


Spalding, Chi, S.W., "The Drag of A Compressible Turbulent Boundary 

Plat® With and Without Heat Transfer", Journal of 
Fluid Mechanics, Voi. 18, Part I, pp. 117-143, Jan. 1964. 


Wallace, J.E. 
Conditions", 
June 1967. 


Turbulent Boundary Layer Studies at Cold Wall 
1967 Heat Transfer and Fluid Mechanics Institute, L** Jolla, CA 


Coefficient Predictions Utll- 

DAC-i9801,^Nov^^ll66^^ * Douglas Aircraft Company, Douglas Report 


Table 5.12 


RAREFIED FLOW SHARP FLAT PLATE 
HEAT TRANSFER EQUATIONS 


‘ ' •< . ' . ?'j 

OF foui QUALIIY 


Rarefied flat plate heating correlations based on the work of Shorenstein 
and Probstein are given below. 


0) 

T, ' 

(2) 

Re« 

(3) 

c* < 

(4) 

$ • 

(5) 

^Hsl 

(6) 

Ch 

4 

(7) 

Heat 




Ww 7 ^ 

Tyj 



(6) l09i«3 + 1.10)j 


for 9 < 0.1 


q • Ch (Ho - 


From Shorenstein and Probstein' 



Nomenclature 

CH Stanton Number 

H Enthalpy 

M Mach Number 

T Temperature 

U Velocity 

X Running Length 

p Density 

y Specific Heat Ratio 
p Viscosity 

Subscripts 

• Free Stream 

w Wall 

^ . 0 Total 

Reference 

Shorenstein* M.L. and Probstein, R.F., "The Hypersonic Leading-Edge 
Problem," AIAA J., Vol. 6, No. 10, October 1968. 
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The rarefied flow heating to sharp and blunt cones based on the work of 
Engel and Praha raj Is as follows: 


(1) T|. • + (Tj ♦ T^^)/2 - Tg cos*6-/3 


(s1n*e + P»cos*e /p^Uj) 


^tY«C★cosec 
(6) Correlation Equation 


logj, (C) - ! (log^gXc) 
1-0 


Sharp 


a# - -0.344074 
a, • -0.349130 
aj • -0.104455 
a, • +0.022766463 


Blunt 


a, - -0.647813 
aj • -0.365587 
aj • -0.0143793 
a, - +0.003281793 


(7) Heat Transfer 

t ■ p. «- Cu (Ho - IL) 




0.9 CH/(slife^ a»*e%/«uUU*) 


Table 5.13 (Cent. 1) 


I s=;;s=Jlli .^igajjigKii ggiB8KgmaBl55ii 


0RiGiw^.L iv.C:^: 

OF POOR QUALITY 


iSrfSiiBsass&gaBaBss 

l!i!i laiiaBsiK 



ST "T ■ mmxm^wua Or^lM^Biaif BUlfea 



— 


*■>■■ MMiaiiiniaiidMraHgnaiaBi 
' • ■ * * Bni Ufa Bu IB Mfi: lua la H ■ ■ ■ 






— 

wsiaiBEipaE^s 

;^:,s|3i:[»iigessssE8ESi^» 

li^sissr ' ■ — 1 1 1 


l^SlgSgSIIB! 


I (»c 

^MiiMijaniB.aeL9ManMRaBiu. K»SMDBaaiM*a 

■ ^ b filsa fif» •; PUiB 1 05 1 aa » ir t ■ 




ttHTSIMitra tt^eiBK 

^tmm3xat%ml ii;- 



s|iu:^ii£^i£gssisggsgsisiigs£:si^i^ 

«{H|MciH^w»Me^p5fS||^SpiiHSSSB8SiBSS&KLjh'^!^8IS! 


laiaBB) 

iNamaMBaBavii 

'^■MlBK8MBilllll 

«!W«BiBiaeiie 


iaiffliasi.. 

mmuK 


t^^ggiyfiseiiEi 

'i£SEgisaginl 




iisipiiniisssoigHti 


iiiiiilili 


Sharp Coni Heat Transfer Data 

Stanton Number 
Enthalpy 

Free Stream Mach Number 

Temperature 

Velocity 

Post Normal Shock Compressibility 
Density 

Specific Heat Ratio 
Viscosity 

Free Stream 
Wall 

Post Normal Shock 
Total 


wIlTECH'iS^'RTRMJ^lt^O.L'^ U78*''"*’ '**"*"**‘ 


0 



Table 5.14 

•p/i." FLAT PLATE HEATING 


ORIGINAL PAGE 19 
OF POOR QUALITY 


The method was developed by Hanks from Boeing and 

by Nagel and Thomas. 


For laminar flow. 


h 


L 


0.332g 


pO.645 V L 


For turbulent flow. 


0.185g 


.0.645 


^ ^r^r 

>'T” [l09,(| (R.r ♦ 


where 


i « 1 (Le°*®^ - 1) 


D,r » f(p^, T^) and L^ « 1.4 


VjM • y, 


/h^Y^^ (T^ ♦ 200) 

yr) 


Nomenclature 

g Gravitational Constant 

h Heat Transfer Coefficient 

H Enthalpy 

L Running Length 

Le Lewis Number 

Kj Laminar Multiplier Factor 

Kj Turbulent Multiplier Factor 

P_ Prandtl Number 

r 


is dc ')tnented 


3000 ))^*®®^ 
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Table 5.14 (Cent. 1) 


Re Reynolds Number 

T Tenq>erature 

V Velocity 

p Density 

M Viscosity 

Subscripts 

e Edge 

D Dissociation 

L Laminar 

r Eckert Reference Condition 

t Turbulent 

T Total 

References 

Nagel et al., "Analysis of Hypersonic Pressure and Heat Transfer Tests 
on Delta Wings With Laminar and Turbulent Boundary Layers," NASA, CR-535, 
August 1966. 

Thomas et al . , "Advanced Re-entry Systems Heat Transfer Manual for 
Hypersonic Flight," Technical Report AFFDL-TR-65-195, October 1966. 
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' - i< 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Table 5.15 

AVERAGE SEPARATED LEESIDE ORBITER HEATING 


The average separated leeside heating relations are based on correlations of 
wind tunnel data for the Space Shuttle Orbiter configuration by Bertin and Goodrich. 
The average leeward heating (turbulent or windward flow) Is given by: 

\ - (1.067^^^+ 0.7905)(0.00282(Re^g)'°-^^) 

^ns * ^"^w^ref^^ns 

Rref • 1*0 foot for full scale 

^ “ pA («t " “w^ \ 

where the windward to total enthalpy ratio replaces the temperature ratio 
In the original paper. 


Nomenclature 


H 


wwd 


'I:*’ ; 
( 

I ‘‘ ■ 

! ‘ 1 

j ’ , 

f ’ 

T; . ■ 1 


«t 

»w 

Re 

R 

P« 

Q 




Windward Wall Enthalpy 

Total Enthalpy 

Local V>^11 Enthalpy 

Reynolds Number 
Radius 

Freestream Density 
Freestream Velocity 

Average Leeside Heating Rate 
Post Normal Shock Viscosity 



l.c«««ard surface area over whicti the h^Mrantftr 
measuremenu were averaged lu obtain 


■ M 

•'1 


t 

. < < 

f 


: S’ 5- 

f ’ I 
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V/V' 




Table 5.15 (Cont. 1) 


OR'niNAL F 13 
OF POOR quality 


A comparison of the correlation and data from the original paper are 
given below. 




A $aumon number for the teeiAtnl **heparmte<l** refloii. 


Reference 

Bertin, J.J. and Goodrich* W.O. » "Effects of Surface Temperature and Reynolds 
Number on Leeward Shuttle Heating," Journal of Spacecraft . Vol. 13, No. 8, 
August 1976, pp. 473-480. 




Oi«GIWA?- F'ACjS Ki 

Table 5.16 OF POOR QUALITY 

PEAK FLAP INTERFERENCE HEATING 


The procedure for calculating peak flap Interference heating follows that of 
Fivel except for the separation length calculations. 



_L 


Flap 

Flowfleld 


TRANSITION 


Re^/MJ < 400 
> 400 


Laminar 

Turbulent 


where 


Undisturbed B.L. edge Reynolds number based on boundary 
layer thickness 

Boundary Layer (B.L.) thickness 
Edge Mach Number 


INCIPIENT SEPARATION (Kessler, Reilly and Mockapetrls) 


Laminar 


Cp^ « Incipient separation pressure coefficient 


’®9loS^ “ '0-361397 log^Q • 0.662427 


Turbulent 


If 2.6 < log|QRe^/MJ <3.8 

Cp^ - 0.2615 (|log,QRe^/M^ - 3.8|)^‘® + 0.405 






o 


o 


Table 5.16 (Cont. 1) 


I, ?-:4! 

OF POOH QUALJTY 


If log^pRe^/M^ >3.8 


Cp^ - 0.0354 Oog^gRe^/M^ - 3.8)’*® ♦ 0.405 


DEFLECTION ANGLE 


p 

2 \ 2-K- (yt) - ly*D^ 

tdn 6|> ) 1 

(y + 1)_^+ Y • 1 



Separation Region Pressure Coefficient 
Deflection Angle of Separated Streamline 


Note If Cp * Cp » Incipient Separation Coefficient 
then 6 q « e^i *> Wedge Angle for Incipient Separation 


PLATEAU PRESSURE (Wuerer and Clayton) 




-1/4 


-1/10 /„2 _ n-1/4 


(Slat), ‘ 


PPLAT 

Re„ 


Plateau Pressure Coefficient 


Reynolds Number Based on Distance along a Streamline 
to the Hinge Line 
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- i 




o 


0 



Table S.ie (Cont. 2 ) 



SEPARATION GEOMETRY 


UKJ'Uw-.v . . V 

OF POOR QUAUi'V' 


i 



Laminar (Giles and Thomas) 


L ,/6 


35 


M 


e L 


(^) 


Re 


,1/4 

ex 


0.98 


Turbulent (Popinskl and Ehrlich) 


where 


Pe 

P 


PL 


X 

6 


T '^ex 




025 




Upstream Edge Pressure 

Plateau Pressure 

Post Shock Flap Pressure 

B.L. Running Length to Hinge Line 

B.L. Thickness Upstream of Separation 


Geometry 

••H 



L^ s1n (6p - 6p)/s1n(180 - 0p) 
Lj s1nep/s1n(180 - 6p) 




•i 





■«« 




u 


o 


Table S.16 (Cont. 3) 


Orra’-"''"’- . ’i 
OF pi’*--!' 


PEAK REAHACHMENT HEATING (Bushnell & Weinstein) 


'shear 


p Uc6 

F s 


M^^s1n(ep - 6p) 


where 


«F ■ 

o*. ■ 
>s. ■ 

<‘.>i ■ 


Velocity at Reattachment Region 
Density at Reattachment Region at Wall 
Viscosity at Reattachment Region at Wall 


6 + 5 


(^) 


(65)^ - 6 ♦ 1.6 Lj /13 


where 


Shear Layer Edge Velocity 
Shear Layer Edge Density 
Shear Layer Edge Viscosity 


PEAK STANTON NUMBER AT REAHACHMENT 


(St)pK - ^st^*^®shear^"^^ 


where 





OP!OilV 

Table 5.16 (Cent. 4) 



^»“3 f. 1 

VM'®,) l‘>'J 


Reattachment for Laminar and 
Turbulent Separated Flows (Bushnell and Weinstein) 
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o.. 


■k 


o 


Table 5.16 (Cont. S) 


REFERENCES 


Five!. H.J.» "Numerical Flow Field Program for Aerodynamic Heating 
Ana»ys1s» AFFDL-TR-79-3128, Vol. 1, December 1979. 

*.!*** ^\^^* Mockapetrls, I.J., "Supersonic 

Turbulent Boundary Layer Interaction With An Expansion Ramp and 
Compression Corner." McDonnell Douglas Astronautics Company, St. 

Louis Report MDC E0264, 17 December 1970. 

Wuerer. J.E. and Clayton, F.I., "Flow Separation In High Speed 
AprflNses ^*'^^**^ Stttte-of-the-Art," Douglas Report SM-46429, 

Giles, H.L. and Thomas, J.W, "Analysis of Hypersonic Pressure and 
Heat Transfer Tests on a Flat Plate With Flap and a Delta Wing With 

Rudders," National Aeronautics and Space 
Administration, Washington. D.C., (NASA CR-536), August 1966. 

Popinski, Z. and Ehrlich, C.F., "Development Design Methods for 

Aerodynamic -Control Characteristics," AFFDL- 
TR-66-85, September 1966. 

Bushnell, D.M. and Weinstein. L.M., "Correlation of Peak Heating for 
Reattachment of Separated Flows . " J^ rna l of Spacecraft and Rockets . 
Vol. 5, No. 9, September 1963. pp. 1111-1112. 
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Table 5.17 

FIN-PLATE PEAK INTERFERENCE HEATING 


Oi^QiNAL FAUE 19 

OF POOR QUALITY 


The peak fin-plate Interference heating can be calculated using the method 
presented by Fivel, which basically came from Hayes. The peak Interference 
heating occurs along a line near the 
fin on the plate. The peak heating 
angle Is correlated using: 


6pk “ 0.24 (6g - Op) + Op 

The location of the peak heating 
line Is given by: 

Ypp ® Xp Tan(6pj^ - Op) 



The boundary layer thickness 
at the leading edge of the 
fin based on the running 
length to the fin leading 
edge, X^, Is required. 



The peak pressure ratio Is based on correlated test data 



(Me S1ne^)”pk 
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0 











Table 5.17 (Cent. 1) 


Likewise, the peak heating ratio Is based on a correlation of test data 



h 


u 


(Me S1n6j -1.0)ng^ + 0.75 


The exponent n and coefficient n^^ are obtained from Interpolation 
of a faired curve through data from Hayes. 


x,/« 



0.0 

1.25 

1.25 

1.50 

2.50 

1.71 

3.75 

1.87 

5.00 

1.98 

6.25 

2.03 

7.50 

2.09 

8.75 

2.14 

10.00 

2.18 

11.25 

2.20 

12.50 

2.23 

13.75 

2.25 

15.00 

2.26 

17.50 

2.29 

20.00 

2.32 

22.50 

2.34 

25.00 

2.36 

27.50 

2.38 

30.00 

2.40 


X 

fie 

o» 

"st 

1.0 

1.62 

2.0 

2.20 

3.0 

2.98 

4.0 

3.57 

5.0 

3.87 

6.0 

4.00 

7.0 

4.05 


If X,/« > 30, • 2.4 

Is obtained by linear 
Interpolation given X./6. 

A 


If Xp/6< 0.183. n^^ • 1.0 

If 0.183 < Xj/6 < 1.0, 

n^^ • 0.35612 ln(Xp/6) + 1.62 

If 1.0 < X./6 < 7.0 

n^^ determined by linear 
Interpolation given X^/6 
If Xp/6 > 7.0 

n^^ ■ 0.35612 ln(Xp/6) ♦ 3.357 


Table 5.17 (Cent. 2) 


Nomenclature 


Subscripts 


Heat Transfer Coefficent 
Pressure 

Edge Mach Number Upstream of Fin 

Distance Aft of Fin Leading Edge In Streamwlse Direction 

Running Length to Start of Fin In Streamwlse Direction 

Running Length Along Fin 

Peak Heating Location Normal to Fin 

Effective Angle of Attack of Fin 

Angular Location of Peak Heating 

Shock Angle 

Upstream Boundary Layer Thickness 


Peak 

Undisturbed 


Hayes» R., '' Prediction Techniques for the Characteristics of Fin 
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Section 6*0 


LB40TB AND TRANSITION ROUTINES 

This seetlon describes the snbrootiaes vhlch deteralne bonndsry layer trss~ 
sltlon, rnanlng length corrections tnd cross flow sdjnstaents* 
d.l SUBROUTINE TRANS AND EDPRAN 

The purpose of routine TRANS is to deterailne If the flow is Isalasr* trea** 
sltlonsl or turbulent* The equations for each of the transition criteria are 
presented in Table 6,1* If a transition option greater than 8.0 is used* then 
an error sieasage is printed from TRANS and the routine returna NERROR * 1 to 
RAIN* This causes MAIN to void the remainder of this problem and read in data 
for the next ease. If the flow is determined to be laminar* then the heat 
transfer coefficient returned from TRANS is equal to the laminar heat transfer 
coefficient. If the flow Is fully turbulent* then the heat transfer eoeffieient 
returned from TRANS is the turbulent heat transfer coefficient* If the flow is 
determined to be transitional* then the heat transfer coefficient returned from 
TRANS is determined by linear i:.terpolation between the laminar and turbulent 
values. 

The interpolation parameter for TRANS “ 1 and 2 is time (TINE)* The inter** 
polatiott parameter for TRANS • 3 is the Reynolds number based on running length 
(corrected for crossflow) or body diameter and free-stream properties. However* 
no correction is made to account for axisymmetrle flow. For TRANS “ 4* the ln~ 
terpolatlon parameter is compressible momentum thickness Reynolds number. 
Transition option 5 uses 

PARA - ^ 

as the iaterpolstion parameter. For transition option 6* the interpolation 


83 


parueur la tU laminar roaning length corrected for croaaflow and corrected 
for aalaymmetric flow (Mangier tranaformatios factor). Transition option 7 nset 
the inpnt running length as the interpolation parameter. Tranaition option 8 
uaea Re ^ /Me as the interpolation parameter. Whenever the flow ia laminar, then 
PCT • 0.0, turbulent PCT • 1.0, and transitional 0.0 < PCT < 1.0. The parameter 
PCI represents how close the flow is to being fully turbulent based on the above 
interpolation parameters. The parameter PCT is also used in the output routine 
(VANOOT) to identify the flow type. It should also be noted that the recovery 
enthalpy for transitional flow is also determined by linear interpolation 
between the laminar and turbulent recovery enthalpies. 


The following is a list of the parameters used in the CALL TRANS statement: 


ENa 

ENCT 

PARAl 

PARA2 

PARA 


TINE 

ENC 

BRE.^V 

NTR 


PCT 


1RFLAG 

EL 

ELFAC 


ENL 

ELL 


(input) L»minar heat transfer coefficient 
(input) Turbulent heat transfer coefficient 
(input) Onset of transitional flow parameter 
(inpnt) Onset of fully turbulent flow parameter 
(output) For transition flag 4, this is equal to the 

compressible momentum thickness Reynolds number. 

For transition flag 5, this is equal to the MDAC- 
EAST transition parameter. For transition flag 6, 
this is equal to the corrected (crossflow and 
axisymmetric flow) running length to the onset 
transition. This is not used for the other 
transition options. 

(inpnt) Trajectory time 

(output) Beat transfer coefficient 

(output) Recovery enthalpy 

(output) N’I'k ** 1 corresponda to laminar flow 

NTR - 2 corresponda to tranaitional flow 
NTR • 3 corresponds to turbulent flow 
(output) Fraction of turbulent flew 

P((r “ 0.0 corresponda to laminar flow 
fCI “ 1.0 corresponda to fully turbulent flow 
0.0 < PCT < 1.0 corresponda to transitional flow 
(Input) Transition option 
(inpnt) Geometric running length 

(input) Multiplication factor which, when multiplied times the 
running length to the onset of transition, gives the 
running length to the onset of fully turbulent flow 
(inpnt) Laminar Mangier transformation factor to account for 
axisymmetric flow corrections 

(input) Laminar running length corrected for crossflow effects 
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Another routine which is closely related to TRANS is EDPARM. Routine ED- 
FARM it need to deter«ine the MDAC~E onset of trsnsition psrsaeter solely ss s 
function of snfle-of-sttsok. Ike sntl«-of-sttsok tenge is froa 0 to 70 degrees. 
There ere sis dets points in the dsts table (o ■ 0» 30. 40* 50* 60* end 70 do** 
gross) along with the corresponding values of the transition paraaeter. 
Interpolation at angles-of-attack not in the table is perforaed with the aid of 
routine IBLIN. 


6.2 SUBROUTINE VRUNL 

Subroutine VRUNL is called by MAIN and is not called unless heat transfer 
technique 3* 4* or 5 (flat-plate laminar and turbulent aethods) is being used. 
Subroutine VRUNL insures that the length parameter used in the turbulent heat 


transfer equations (when using transition options 3* 4* 5, 6* or 7) is corrected 
for the assumption that the running length should be measured from the onset of 


transition. The ergument of subroutine VRUNL contains the following parameters* 


TRFLAG - (input) Nay have values 3.0* 4.0* 5.0* 6.0* or 7.0 

corresponding to transition options 3.0* 4.0* 5.0* 
6.0* or 7.0. (It should be noted that VRUNL will 
not be called unless a positive value is stored in 
location 29.) 

ELIRAN - (output) Running length at onset of transition with or 

without crossflow correction 

Bi.f . - (input) Equivelent running length with or without crossflow 

correction for laminar flow 

ELT - (input) Equivalent running length with or without 

and crossflow for turbulent flow and is corrected in 
(output) VRUNL to account for virtual origin from the onset of 
transition 

ELTP - (output) Equal to input value of ELT 

PARAl - (input) This parameter denotes the onset of fully 

transitional flow 

PARA2 - (input) This parameter denotes the onset of fully turbulent 

flow 

ENL - (input) Laminar multiplication factor 


If PARA2 < PARAl* then PARA - PARA2. If PARAl > PARAl* then PARA - PARAl. 
If IIFLAO < 3.0* then MAIN does not call VRUNL. Therefore* VRUNL is really only 
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B 0 snln|fal for TOFLAQ - 3.0, 4.0, 5.0, 6.0, or 7.0. The valooi TRFLAC “ 1.0 or 
2.0 corrosposd to tino dopendont transition and ara handled tbroogh the main 
rostina. If TRFLAG - 6.0, than VRUNL oalU EDPARM aa a fnnetion of 
angla-of-attaok to dataraina PARA. Tba valoa of ELTRAN is given by 

ELTTIAN - 

If IRFLAG • 7.0, than the logarithm (base 10) of the transition Reynolds number 

RELG, based on edge conditions, is dataralnad as a fnnetion of edge Mach number 

H from a ouilt-in table (North American Rockwall technique). The value of EL** 
0 

TRAN is given by 

ELTRAN • (10 V*e^/Pe'‘c 

The interpolation to determine logj^^Re^ is performed with the aid of routine 
TBLIN. If TRFLAG • 5, then ELTRAN is the same as for option 5. However, the 
value of PARA • PARAl is input through the routine argument to define transition 
onset instead of being determined by EDPARM. If TRFLAG > 4, then the running 
length at the onset of transition is determined based on the compressible momen- 
tum thickness Reynolds number: 


ELTRAN 



.664 


1 


a a 

p p u 

e 


where PARA ■ Re* is input to define the onset of transition. If TRFLAG ® 3, 

6 

then the running length at the onset of transition is determined bssed on the 
Reynolds number computed using boundary-layer edge conditions: 


where Rej^ . para. 


ELTRAN 


RClP 


P.u 
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After the running length at the onaet of traneition ia calculated, the tur- 
bulent running length ia atored (ELTP - ELT) and a check 4a Bade to detersine if 
ELL < ELTRAN. If ELL < ELTRAN, then ELT • ELTP. If ELL > ELTRAN, then ELT • 
ELT - ELTRAN. Thue when VRDNL returna to the call log routine, ELTP containe the 
input value of El.T and ELT may or nay not contain the input value of ELT depend- 
ing on whether ELL < ELIRAN or ELL > ELTOAN. If BLL > ELTRAN, then ELT conteina 
the turbulent running length seaeused fron the onset of transition. 
d.3 SUBROUTINE CSSFLW 

Subroutine CRSFLW is called by MAIN and the purpose of this routine le to 
account for the effects of stresBline divergence on basic geoBCtrlea at 
angle-of-attack. IVro types of geometry are considered in this routine: a con- 

stant width shape (rectangle) and a sharp-edge triangle (delta wing). For each 
of the two geometries it is possible to get an ideal or reel gas calculation for 
both ItcJnrr and turbulent flow. In each circumstance, a correction factor is 
applied to the input geomtric length, and this corrected running length ia used 
in the flat-plate heat transfer e<,uationa in place of the input geometric 
length. The parameters in the argument of the CALL CRSFLW statement are lieted 

and defined below: 

CFFLG - (input) T^pe of crossflow correction selected 

• Constant width rectangle assuming ideal gae 
chordwise velocity gradient 

• Constant width rectungle using a real gas 
chordwise velocity gradient 

* Sharp-edged delta configuration assuming ideal 
gas chordwise velocity gradient 

* Sharp-edge delta configuration using a real gas 
chordwise velocity gradient 

ELL - (output) Laminar running length corrected for crossflow 

ELT - (output) Turbulent running length corrected for 

EL - (input) Ph /ical geometric running length input into MINIVER 

ELMBDA - (input) Delta sweep angle 

DSDBO - (input) Rectangle width 

CORNR - Unput) Rectangle corner radius 

BE - (input) Velocity at the edge of the boundary layer 

PE - (input) Pressure at the edge of the boundary layer 
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/iLPHA 

DDOT 


( input ) 
(input) 

( input ) 
(output* 


- (input) 


XMAtaiU - (input) 


Pontity nt tho odge of the boundary layer 

Presaure upstrean of the previous shock (if no shock, 

then freeatream pressure) 

Surface effective angle-^of-attack 
Designates the nondinensional crossflow, 
stagnation-point, velocity gradient at a point on 
the centerline (of the wing) as a ratio to the 
velocity gradient at the stagnation point of a 
sphere with a diameter equal to the planform width 
Velocity upstrean of the previous shock (if no shock 
then free-strean Mach number) 

Mach number upstream of the previous shock (if no shock 
then freestream Mach number) 


Having entered GtSFLW, the first check that occurs is to determine if ALPHA 
1 0, if so then ELL * EL, If ALPHA > 0, then the decision is made as to which 

crossflow option should be used. If CFFLG - 1, the necessary program input is 

CORNR^and DSUBO. The routine then proceeds to compute the parameter ^ 

TP ^ ^ crossflow is used. This indicates the velocity gra- 

dient is so small that crossflow is insignificant. If ^>6 . then the 

Oq 

routine shifts to a simplifiet’ calculation of ELL and ELT (this situation cor- 
responding to the case of a very large velocity gradient). 

TV -R XV i , 

10 ° < vj— < 6*“®® ^*^® ctendard crossflow correction expressions *ie u-ved 

to compote ELL and ELT and then CRSFLW returns to MAIN. If CFFLG » 3, then the 

ideal gas crossflow correction for a sharp-edged delta is used to compute ELL 

and ELT. When CFFLG = 3, the only necessary Input is the delta sweep angle, 
ELNBDA. 


If CFFLG » 2, then the crossflow correction for a constant width rectangle 
using a real gas chordwlse velocity gradient is used to compote ELL and ELT, 
The required input are DSUBO and DDOT. If CFFLG = 4, then the crossflow correc- 
tion for a sharp delta configuration using a real gas chordwise velocity gra- 
dient is employed to compute ELL and ELT the required inputs arc the delta eep 






nor«»l conpoMot of opstrean Maoh aombar XMACHD ia leaa than ona. If XMACBO < 
1* than a anbaonlc tachnlqna ia naad to conpota tha ocoaaflow ooxraetion. If 
XMACBD 1 1, than a anpaxaonlo taobnlqoa ia naad for oroaaflow eorraotion. 
Howavax for CFFL6 ■ 1 or 3, no diatinotlon ia sada for XMACBU < 1 and 
XMACBD i 1. 
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Table 6.1 

TRANSITION OPTIONS 


Of the^followiiig^***^*"^^^ transition options In the LANKIN code consisting 


1. Time Dependence: Laminar to Turbulent 


t < tj 
t > t 


II 


Laminar 

Turbulent 


ti < t < tjj Transitional 

t-t| Percentage of Fully Turbulent Flow 

^ ° ‘T ' . ' t ' 

^II ^I 

2* Time Dependence; Turbulent to Laminar 


t < t 


I 


Turbulent 
t > tjj Laminar 
ti < t < tjj Transitional 

Percentage of Fully Turbulent Flow 


*irh 


3. Reynolds Number Dependence 

*^®oo 0 “ fo<* Swept Cylinders 

for Plate Options 




Re < Re 
Re > Re 


I 

II 


Laminar 

Turbulent 


Rej < Re < Rejj Transitional 

Re-Rej Percentage of Fully Turbulent Flow 
• TliJ-pej 

Compressible Momentum Reynolds Number Dependency 

0 0.664 Tp y Re^ Compressible Momentum Thickness 

\ o u T Reynol ds Number 
f ^"e*^e L 


on 


Table 6.1 (Cont. 1) 


R6q < Regj Laminar 

* 

Reg > R®0|| Turbulent 

Regj < Re* < Regjj Transitional 

* 

^ _ ROg - ROgj Percentage of 
^ ■ W, 


'01 1 


■Re 


ei 


Fully Turbulent Flow 


5. MDAC-E Transition Parameter 

This transition option was developed by Nasek and Kipp 

ip • Re* / [Me(pgUg/Ug)“] 


♦ 5 ’f'l 
^ > ’I'll 

\lil < il> < 4>jj 


n 




Laminar 

Turbulent 

Transitional 

Percentage of Fully Turbulent Flow 


6. MDAC^c. Transition Table Look-up Dependency 

The built In tabular function of angle of attack from Masek and Kipp 
was developed for conical flews. Transition Is based on the parameter 
^j/ from option 5. 


a 

^O9io ’I'l 

0 

1.0 

30 

1.0 

40 

1.0412 

50 

1.1363 

60 

1.30103 

70 

1.69897 


Linear Interpolation 
Is used 


Once the transition onset parameter. Is calculated, the transition 
length Is determined: 



.6 
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Table 6.1 (Cent. 2) 


If a crossflow option has been flaqqed, the followlnq procedure Is used: 
Since the transition parameter ij' was derived for complex shapes usinq 
the laminar crossflow length In the calculation of Re*, It Is necessary 
to remove this correction factor so that the resulting transition length 
can be compared with the geometric length. 


For crossflow options for rectangular shapes the transition length is 



-0.5 V In 



I 


For delta wings options the transition length Is 


•■TR “ ®4r' 

When the Mangier transformation, Tj^, Is used, the transition length Is 


•■TR “ 4R' 

Two methods for treating the extent of transition have been provided. 
The first Is rather simple in that the fully turbulent length, Lpj, Is 
a factor times the onset length ^ 

•■FT ' '^TE ’ '■TR 


The second method defines the factor In terms of the boundary layer 
edge Reynolds number at the onset length, which Is defined as 


“e^e^TR 


RRte ' 


The table, which Is built In, consists of 


iog,o ««te 


s. 

5.64 

6.518 

2. 

6.778 

1.702 

7. 

1.605 

7.301 

l.n35 

10. 

1.535 


( Linear Interpolation 
Is used 
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Table 6.1 (Cent. 3) 


The state of flow Is detentilned by 


‘■LCF 

- ‘'TR 

Laminar 

kcF 

- ‘■FT 

Turbulent 

••TR 

‘•LCF "" ‘■FT 

Transitional 

n ■ 

‘•FT “ ‘■LCF 
‘■FT ’ ‘•TR 

Percentage of Fully Turbulent Flow 


Since the transition criteria was developed for e conical shock flowfleld, 
some Judgement must be used In applying It to other situations. 

7. NAR Transition Parameter 

A criteria developed by North American Rockwell Is based on edge Reynolds 
number as a function of edge Mach number. The built In table Is 

^ ^O9l0 ^^TR 

0.5 5.30103 \ 

1.0 5.54407 I 

ko 6*00^^^ i Interpolation Is used 

3.5 6*. 00 ! 

4.5 6.0791P I 

5.5 0.25527 I 


• Re.|.|^ Transition length 

The transition extent Is the same as for option 6. The state of flow 
det'-rmlnatlci. Is based on rather than Lj^gp. 

8. Compressible Reg/Me Dependency 


Table 6.1 (Cent. 4) 




• Re^ /Me 

Rf < 

Rfl 


Lami nar 

Al 

OC 

Rfii 


Turbulent 

Kfi < 

:Rf < 

Rfl I 

Transitional 

n « 

iL 

■ Rfl 

Percentage of 


Nominal Values 
Rfl ■ 150 


Rfll - Rfi 


Nomenclature 




S 

Cross Flow Parameter 

m 

V 

Chordwise Velocity Gradient 

1 

Length 

X 

Cross Flow Corrected 

M 

Mach Number 


Running Length 

t 

Time 

p 

Density 

U 

Velocity 

y 

Viscosity 



T 

Mangier Factor 

SuDicrlpts 




0 

Diameter 

L 

Laminar 

e 

Edge 

LCF 

Laminar Cross Flow 

Geo 

6ec*ietr1c 

TE 

Transition Extent 

I 

Initial 

TR 

Transition 

II 

Final 

e 

Momentum Thickness 
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Section 7.0 

APniCATION TECHNIftOES 


n. fuioM .f tM. 1 * «» •' “• •**“**' “ 

•ppirUi tk. k"*™ *• “• “ '“**‘'‘* *■ ** 
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LAmilk. B...< .. tk... ....It. .1 P*.t ..P..1.... .Itk tk. . ..t »I 

cooacndcd options nco giyon. 


7.1 GENERAL APfLICATION 

In order to nchieve the beet prediction nethod nelng the LANMIN code eever- 

nl key pnremetere need to be examined. 

1. D.t.nil>. tk. A..k .ptlo. .kick ..t.kll.k.. tk. kotmtorp Up.. 

t Djt*»"M*S.*iIoiit.l« N..|l.. tl...fCM.tlo. .kick jl... tk. 

* nnnroprinte rimning length ndjnatnent for the given geometry, 
s. mJSii: thn**.!~* option which beet metche. the enrfece 

4. tk. k..tl.» ..tko4 ..4 

factor which beat fits the experimental heating data. 


Ik. .kck optic. 1. MMllp MUctcd kp 4.M»l.l.f .kick cptlo. p.c4.c.c 
tkc .ppcc.l>.t. rtcck ...ic .kick P.OCC..C. tk. ko..4.«r Up.. ««. nl. 1. 
nsnally eatabliahed by using wind tunnel achlieren. 


Ik. Il..tl.. t....f 0 M.tl 0 . ...4 1. 1. tk. t. 0 Mt.t. t....fc...tlc. 

c( tk. .U.1.1 u.ttk lo. . ko4p l.to tk. ...lo.U.t Ic.tk fo. • fUt pl.t.. 
Ik. ..tk«.tlc.l ..Utlc. .kick ... k. ...4 . c<»po«.4 ,.c.t.p k.4p .kick 


rs 


t 

r s 

B ^ ft 

1 ' 

^ 1 

’ j* pd» 


f r*ds 


Laminar 

9S 

Turbulent 

) 

-* 




:r 





where r le the redine of the body norael to the exit of eyanetry end e i) the 
enrfeoe dletanee from the noee to the point of internet on the body* The 
preeedini equations can be nnaerieelly Integrated for compound body shepea to 
obtain the laminar and turbulent Mangier factor aa a function of running length. 

Selection of the preaanre option and Reynolds analogy factor for a problem 
ia baaed primarily on experience* In order to document some of this experience* 
several sample cases will be examined* Figure 7*1 shows a comparison of wind 
tunnel pressure data with the LANHIN tangent cone option and a method of cha«'ae*- 
teristic (MOC) caleulation* The data and theories are for the cylinder section 
of the Space Shuttle external tank at xero angle of attach* For this situation* 
the tangent cone option does not predict the data well* Thia option would also 
correspond to the modified Newtonian pressure option* The MOC results do com* 
pare well with data* Thus* results from a more exact theoretical method should 
be used as input to LANMIN, where available* to establish the edge pressure* 

Beating results from LANMIN using both pressure options are sho*)fl in Figure 
7*2 using the best heating prediction method* The overprediet ion trend of the 
tangent cone option* shown in Figure 7*la* is translated into an over prediction 
in heating as shown in Figure 7*2* Likewise* when the pressures from the more 
accurate MOC solution are used* the heating is more accurately predicted* 

The Reynolds analogy factor effect is examined in Figure 7*3. The Colburn 
aualogy factor ia shown to overpredict the data and a unity Reynolds analogy 
factor underpredicts the data* The von larman analogy factor when used with the 
Spalding‘>Chi heating sisthod with accurate pressures provides the best comperison 
with data. Thia trend le substantieted further by the results given in Figure 
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7. 4, Tb* prettnre at tmill to aoderate anglot of attack aloog tko wladward 
ataaaalina ia pradictad wall by both tha tangaat aona option and HOC aolntiona. 
Thna tha haatinp data in Fifura 7.4a aoaparaa wall with ealanlationa nalnt both 
praaanra optiona. lha theory and data aoapariaona in Figure 7.4b and a doaon* 
atrata that tha won Earaan Raynolda analogy factor givaa the beat pradietion. 

Saaulta froa tha Spalding-Chi option naing tha yon Karaan Raynolda analogy 
are ahown in Figure 7.5 along with an independent heating theory for a different 
geoaetry. Heating data for a aphare~aone**cylinder (SRB) body are ahown. The 
DIRLIN reaulte are froa an integral aoaentua solution baaed on the work of Dir~ 
ling (Ref. 7). The two theoriea are in good agreaaent and eoapare wall with 
the wind tunnel data. 

Another application of the LANMIN code ia ahown in Figure 7.6. Thia caae 
ia a eharp wedge at three anglaa of attaek. Thraa-diaanalonal effecta on the 
wedge produced preaaurea higher than the two-diaenaional wedge preaaure option. 
By naing the wedge ahoak entropy and inputting the aeatured preaeuree. the LAN~ 
reaulte eoapare well with data and reaulte froa an independent theory of 
Dlrling (Ref. 7). 

All preceding reaulte and aoapariaona haye been for turbo? ent flow. The 
laainar aethod of Eckert ia LANMIN ia eoapared with wind tunnel data i*' Figure 
7.7. The tangent cone preaauze ox>tioa waa need for thia caleulation. Being 
aore preeiae preaaurea would haye alightly lowered the predicted yaloea in the 
0.08 ( X/L ( 0.22 and alightly increaaad the yaluea in the range 0.22 4 
X/L < O.dO. 
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Other coaperlsont of date and theory have been preeeated la the tablea of 
Seetloa 5*0» Theae reaalta aloat with other esperleaee were aeed to deteralae a 
reooa»eaded aet of heatiaf optioae* These optioae are givea la Table 7.1. 

7.2 GENERIC ORBOER APHICATIONS 

Beatiag to a geaerio arbiter geometry eaa be ealealated by a aaoeroae coai'* 
biaatloa of optloas for a speelflo loeatloa oa the vehicle. The soeeese la ac** 
®®**l*ly predictlag the heatlag to aay givea looatloa Is deteralaed la large 
part by the proper seleetloa of optioae aad the eagineerlag eklll of the aeer. 
A geaeral approach la givea herela for orbiter appllcatloa. Specific geoaetries 
aad flow eitaatioas may reqaire Bodificatloa of the gaidellaes givea. 

The approach takea follows the eelectioa of the four key paraaetera givea 
la Seetloa 7.1 as apeeifieally applied to a generic orbiter shape. Pigare 7.8 
gives the selected geaeric orbiter deeiga aethods ia eoaaary fora. TTie approach 
takea was to select aethods which woald beet represeat aeaeared data. If a 
deeiga coaservatisa is to be iatrodaced. it ahoald be iatrodaced as a kaowa conn 
staat aaltiplier. Ia order to deaoastrate aad help ezplaia the options eflect" 
ed. a discasaioa is givea of three of the priacipal arses oa aa orbiter. 

Nose Reaioa 

The heatlag to the bottom ceaterliae of am orbiter aoso eaa be ealealated 
®*i®f ihe optioae ia LANHIN by deteraiaiag aa effective sphere valae. The ef- 
fective sphere valae is ealealated with aorael shock eatropy. The heatlag 
aethode of Lees aad Detra— Hidalgo are aeed for laaiaar aad tarbaleat flow 
raapeotively. Coaalder the laaiaar case as aa ezaaple. The heat traaafer eoef- 
ficieat ia deteraiaed by 



( 


ho it tht P»f *Bd Riddoll otlM for radios Ro 
r it tht lootl radios of tho body 

h/ho it tht Lttt dittribotioo ratio at tha local body angl* aosla 
of attach 

An axanpla of this aathod being applied to a tphara-'cont gaoaiatry ia gioan in 
Fig. 7.9. Tha Laaa diatribotion Taloe on tha cone ia avaloatad at 
90 - 20 - IS ■ S5 degrees, and the loeal radios, r. is a fosetion of axial loea- 
tioo. This rather tinple method provides ooite good agreenent with the data 
shown. 

This aathod was also applied to the note region of an orbiter eonf igoration 
a* shown in Fig. 7.10. The area of application was 0 X A/L < O.OS on the or** 
biter bottom centerline where theory and data agree for all three angles of at** 
tach. 


Thos by inpotting the body angle and loeal radios an effectiwa sphere dit> 
tribotioB waloe is ealcolated which repreaents the note region bottom center* 
line. The effect of angle of attack ia aecoonted for by the program caleolating 
an effeetiva body angle. 

Bottom CgtttffUM BfAtlM 

In order to select a osable method for LANMIN for the bottom centerline of 
an orbiter. each component of the ealenlation most be considered. The compari* 
sons given in Fig. 7.10 for X/L > O.OS are sample trials of different options. 
At oat angle of attack, a given set of options may agree with data bnt not agree 
with data at other anglaa of attack. This approach eonfonndt tha affects of tha 
different options in a set. Thna. final agrasaont for one case docs not assure 
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•pplieabilltj ov*r a wida raafo of oonditiona. 


Lot os coatldtr tlia shock shape which controls the boundary layer edge con- 
ditions. Pi|ore 7.11 presents a correlation of shock shapes first published for • 
delta wing data by Dunavsnt (Bef. 14). The generic orbiter data from four 
sources vere added in the present study. Virtually all of the data agrees with 
tangent cone shock theory for angles of attack fron 10 to 45 degrees. Most all 
orbiters fly in this angle of attack range during the highest heating period. 
Thus the shock option to be used is the tangent cone option. 

Determining the correct running length to use is the most complicated and 
most critical decision to be made. The LANMIN code has cross-flow correction 
options and Mangier transformations to assist in correcting the actual running 
length to an eqniwalent running length.* For some geometries these running 
length corrections can be quite accurate and easy to select. This is not the 
case for generic orbiter geometries. 

Another approach ia selected here based on the modification of e».rlier work 
by Pond et al (Bef. 15). iminar heat transfer data are used to calculate -the 
effeetiwe running length to the body location of interest. This procedure finds 
the actual flat plate equivalent running length which is required for input into 
LANMIN. This effective running length is then used for flight prediction calcu- 
lations. A sample case is described to provide details of the approach. 


* Other analyses have developed complex streamline tracing methods which re- 
qaires an accurate and smooth pressure distribution input. These methods are 
out beyond the scope of the LANMIN approach. 
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Pond «t *1 (R»f, iS) t 0 $fd th* »bown in Pif. 7»l2 to obtain 

bontiag data. Tba boating data war# naad with Bokarta boating natbod to aalat*** 
lata tba a^iniyalant running langtba aboun in Fig. 7.13. For tbia gaonatrj it 
ia auidant that tba aquivalant running langtba ara a strong function of angla of 
attaek and axial location. 

Work by Dunavant (Raf. 14) wbieb vas latar aaplifiad by Nawnann and Ban- 
froa (Raf. 16) daalt with obtaining an outflow corraction for tba running 
langtb. For bypartonie flow, tbia correction can sinply be written aa 

Xa ■ X(tan s/tan a) (7.1) 

Xa ■ aquiwalant running langtb 
X “ actual running langtb 
a " angla of attack 
a • 90 sweep angle 

By looking at tba gaonatry in Pig. 7.12 a single walua to use for c is not awi** 
dent. Howsyar* by taking tba data in Fig. 7.13 and using it in tba preceding 
agination* s can be calculated. The raaults ara abown ia Pig. 7.14. All of tba 
20. 30 and 40 dagraa angla of attac«. data collapse to giya an affactiya sweep 
angla eonplanant aa a function of axial location. Tba a * 10 data ia for aa inr 
flow condition and ia thus not axpactad to correlate. Ike data for a 55 asbi** 
bit a trend wbieb ia not conplataly understood and nay be due to lack of quality 
paint data for tba original beating data. Since tba preceding relation dona 
collapse tba data in tba prinary range of intsrast. it can be used along with 
United wind tunnel data to dayalop aquiyalant running lengths for input. 

To conplats tba explanation of tba procedure, tba aquiyalant running langtb 
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data is Fig. 7.13 was naad in LANNIN for a ooaplataly diffarant orbitar abapa. 
Tha raaalta froa LANMIN ara aoaparad io Fig. 7.13 with data froa a phaaa B or- 
bitar gaoaatry. Tba coaparlaon la qoita good and significantly iaprovad ovar 
tba coaparisons aada is Fig. 7.10. Nota that t!»a tangant wadga praaanra option 
was naad. Unf ortunataly. pratsora data ara not availabla for coaparlaon. 
Bowavar, Mach nnabcr data on orbitar abapaa ara naoally ovarpradictad whan wadga 
pxataores ara naad. Thus, wadga prasauraa ara most libaly too low. This might 
azplain why tha thaoratical resnlts ara slightly lower than measnrad in Fig. 
7.15. 


Thns far only laminar aqnivalant lengths and rasolts have bean discnssed. 
Tha laminar hasting rates ara more sansitlwa to tha aqnivalant length 
(h - Xa*^*) than turbnlast hasting rates (h ->» Xa**). The laminar and tnrbnlant 
equivalent lengths ara diffarant. For example on a cone tha diffaraneas ara: 

Xa, • X/3 

** from Mangier relations 

Xa^ - X/2 

Thus an adjustment must bo made to tha laminar equivalent length to obtain 
a turbulent equivalent length. Since no data ware available for turbulent 
equivalent lengths, the cone ratio was used 

Xo^ - 3Xe^/2 

to obtain input to LANMIN for the ease shown in Fig. 7*16. Tha comparison of 
data and theory is quite good. Note that a 50% increase in the running length 
when used in turbulent heating equations only reduces the heating by 8 percent. 

Tha procedure to use tha equivalent running length approach is as follows: 

(1) Obtain laminar wind tunnel data for the generic orbitar of interest. 

(2) Use LANMIN parametrically to generate laminar equivalent running lengths. 
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(S) t£ dat* ara availa^la, rapaat atap (1) and (2) to rotltj tka laninar to 
torbnlaat tranaforaatlon. 

(4) If lialtad data ara aaailabla* naa atoatioa 7*1 along nltk an affaetlaa a 
baaad on tka lialtad data to aspand tka angla of attaak ranga* 

(5) If no wind tnanal data ana aaailabla naa Pig. 7.13 os eloaaat gaoaetsy 
aqnivalant iangtk data (Nota a aat of plota/tablaa of aqniaalant langtka 
for a aat of gaoaatsiaa akonld ba dasalopad aa a data baaa). 

(6) Uaa tka a^niaalant snnning Iangtk in LANMIN for fllgkt applicationa along 
nitk tka aaaia akook. psaaansa and kaating rata options noad to ganarata tka 
aqnisalant snnning langtka, 

Tkis natkod ka^i baaa disenaaad in tasas of tka botton cantasliaa. Bowavar. it 
kaa a anak bsoadas appliaation. Tkla aatkod aay ba appliad to any body loaatioa 
vkiak ia net iaflnanaad by sapasation os tkoak intasfasanaa. 

2ipg Laadina £dli 

Tka stagnation liaa of a wing ia appsoxiaatad nail by anapt aylindas tkao'* 
sy. Tka skoak liaa pasallal to tka ning onas a aonaidasabla portion of tka wing 
aa akown ia tka sakliasan of Pig* 7.17. Tka Isainas flow data agsaa nail nitk 
tka aossalatioa ia LANMIN. A aoapasiaon for tnsbnlaat flon of 
Baafcnitk-Oallagkas nitk data is ginan in Table S.da. 
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Table 7.1 


RECOMMENDED METHODS 


(iEOMETRY 


HEATING OPTION 


HEAT TRANSFER 
OPTION NUMBER 
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CONE, WEDGE, 

OGIVE & ORBITER 
BOnOM CENTERLINE 
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Fay and Riddell 


Laminar: Lees 

Turbulent: Detra and Hidalgo 


Laminar: Correlation 

Turbulent: Beckwith-Gallagher 


Laminar: Eckert 

Turbulent: Spalding-Chl with 
Von Karman Reynolds Analogy 


Bertin and Goodrich 


Bushnell and Weinstein 


FIN INTERFERENCE 
ON PLATE 


Fivel 
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Fig. 7.1 
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Fig, 7.6; Heating Rate Versus Length for Three Test Conditions at the 
Ames TPS Ca^'^ratlon Test 
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Fig. 7.8 LANMIN Generic Orbiter Design Methods 
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Fig. 7.9 Windward Centerline Heat Transfer Over a 20-0egree 
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Fig. 7.11 Shock-Wave-Angle Correlation 
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Fig. 7.15 LANNIN Theory Compared With Phase 6 Orbiter 
Configuration Bottom Centerline Laminar Data 
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Fig. 7.16 LANMIN Theory Compared With Phase 6 Orbiter 
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